SCATTERING POLES NEAR THE REAL AXIS FOR TWO STRICTLY 

CONVEX OBSTACLES. 



ALEXEI lANTCHENKO 



Abstract. To study the location of poles for the acoustic scattering matrix for two strictly 
convex obstacles with smooth boundaries, one uses an approximation of the quantized billiard 
operator M along the trapped ray between the two obstacles. Using this method Gerard (cf. 
[9] ) obtained complete asymptotic expansions for the poles in a strip Im z < c as Re z tends 
to infinity. He established the existence of parallel rows of poles close to ^ + ijS, fc S Z, 
j G . Assuming that the boundaries are analytic and the eigenvalues of Poincare map are 
non-resonant we use the Birkhoff normal form for M to improve his result and to get the 
complete asymptotic expansions for the poles in any logarithmic neighborhood of real axis. 



1. Some known results on the localization of resonances. 

1.1. Introduction. Following [26] we denote by — Ar^+i^q the Dirichlet Laplacian on a con- 
nected exterior domain M"+^ \ O, where O is compact with a C°° boundary. Then the resolvent 

Ra{X) = (-Ar„+i\o - A^)"' : ^^(Rn+i \ ^ h'{W'+^ \ O) H^iW'+^ \ O), ImA < 0, 

continues meromorphically across continuous spectrum ImA = 0, to an operator 

i?o(A) : i?omp(K"+^ \ O) ^ HUR-+^ \ O) n HI,,M"+' \ O). 

Here H^{W^+^\0) is the standard Sobolev space, Hl{«^+^\0) is the closure of C;;^^p(R"+i \0) 
in the iJ-'^— norm and by L^o^^p, where by .Z^comp mean the elements from that are zero 
outside some bounded set, and by H^^^, H^^^^ functions that are locally in these spaces. 

We recall that RoW is globally meromorphic in A e C when n + 1 is odd and in A £ A, the 
logarithmic covering of the complex plane, when ri + 1 is even. 

Definition 1. The poles of Ro are called resonances or scattering poles. 

Note that we adopted a convention that the resonances Xj are in the upper half plane Im A > 0. 

We cite some results following the review by Zworski [5S] (note his convention that Im A^ < 0.) 

An obstacle O contained in a ball B{0, R) is called non-trapping if there exists T > such 
that every broken characteristic ray (see e.g. [11], Section 24.2) starting in (R"+^ \0) O B{0, R) 
leaves that set within time T. 

Propagation of singularities for boundary value problems, established by Melrose-Sjostrand 
and Ivrii (see Chapter 24 of [TT] and the references given), gives the basic results: 

O is non-trapping, dO is C°° ^ 

(1) there exist a,b > such that a region {z; Imz < alog(l + |z|) + b} 
contains no poles; 

(2) VA^, tt{Aj- : llmAjl < Ariog(l + \Xj\)} < oo. 
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The situation is somewhat more satisfying in the analytic case where one can apply the prop- 
agation of Gevrey 3 singularities due to Lebeau [TB]. In fact, Bardos-Lebeau- Ranch f7 (and also 
independently Popov 19J) used the work of Lebeau to show that 

O is non-trapping, dO is real analytic ^ 

(3) there exist a,b > such that a region {z; Imz < a(l + |z|)^^'^ + b} 
contains no poles. 

Then clearly, 

(4) 3B,C>0: tJ{Aj : |Im | < B|Re A^f < C. 

This is known to be optimal in the following sense [2]: a non-degenerate isolated simple 
geodesic 7 of length on the boundary of an analytic strictly convex obstacle in n + 1 

odd, generates infinitely many poles in any region 

{z: Imz<-P-|z|U, S > tj-, = 2"^Ci cos f^V / p^^^(s)ds, 

\6/ Jo 

where —(1 is the first zero of the Airy function, p is the curvature of 7 in R"+^ and s the arc 
length parameter. 

For strictly convex obstacles the results on pole counting are expected to be more precise and 
in particular the density of poles near the real axis has already been estimated. To motivate 
them let us consider (in [5^ the convention was that Im A^ < 0) 

iVg(r) := tt{Aj : A^ is a pole of i?o(A), \Xj\ < r, < arg Aj < 6*, n - < arg Aj < tt -|- 0}. 

Then for strictly convex smooth obstacles we get 

(5) Ng{r) ^ O{0^/^y'+\ r > 1. 

By applying complex scaling " all the way to the boundary" which is possible because of strict 
convexity, but still using a form of functional calculus from ^U\, (O was obtained in [3T] for 
an arbitrary convex obstacle provided r > r{9). It was noted by Harge and Lebeau [7| that a 
particular choice of the angle of scaling yields, among other things, a cubic pole free region of 
the form ^ for strictly convex obstacles with smooth boundaries. It had been widely believed 
in applied mathematics and known already in dimensions two and three (see [T]). 

1.2. Trapping obstacles: results of Ikawa and Gerard. We consider scattering by two 
strictly convex obstacles. Note that scatterer is always trapping if it is not connected. 
We suppose that 

C = r^i U C M"+\ with C°° boundary dfl,, i = 1,2, 

(6) where fli, fl2 are compact and strictly convex, fii n = 0. 

Let d be the distance between fli and O2 and g dfli the points on the boundary such that 
\ai — = d. Then there is one trapped broken characteristic ray connecting oi and 02. 

Under these assumptions on fl Bardos, Guillot and Ralston (cf. [3], for n -|- 1 odd) show the 
existence of an infinite number of resonances in {z; Imz < elog |z|} for any e > 0. 

Thus their result shows a difference in location of resonances between cases of trapping obsta- 
cles and of non-trapping obstacles. 

The most complete results on location of poles were given by Ikawa and Gerard. There results 
shows that uniformly in the strip < Im z < c, as Re z goes to infinity the resonances are well 
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approximated by pseudo-poles of [3]: 

n 

(7) A„,fe-fc^ + ^^lnz.,.(2a, + l), (fc,a)eZxN", 

i=i 

where Vj are the eigenvalues of the linear Poincare map DK(ai,0). 

Though Ikawa and Gerard only consider the case of odd dimensions 7i + 1, their results are 
valid in even dimensions n + 1 > 2 as well. 

In [T3] Ikawa obtained the first string of resonances closest to the real axis (a = 0) and in (cf. 
[S]) Gerard got the complete asymptotics of all strings. 

We start by defining the canonical transformation of the billiard k, : T*dili t-^ T*dili following 
Gerard For p £ T*dfti in a neighborhood of (ai,0) G T*dfli we draw an outgoing half-ray 
7 issued from p and we denote 7' the half-ray reflected from dfl2- If 7' intersects dfli then we 
define k{p) the projection on T*di^i of the point of intersection between 7' and 9T*(R"+^ \ ili). 

For p S T*dn,i, close to (ai,0), we denote by Ki{p) the projection on T*dVl2 of the point of 
intersection of 7 and 9T*(R"+^ \ ^12)- ki is the canonical transformation from T*dili to T*dVl2- 
If we define in the similar way the canonical transformation K2 '■ T*dfl2 1— > T*dfli then the 
billiard map is the canonical transformation 

(8) K = K2 ° f^l- 

According to Petkov ([U]), [H], Corollary 2.3.3) and Bardos, Guillot, Ralston ([5], Proposition 
3) using that fii, rt2 are strictly convex, the eigenvalues of Dk(/9i,0) are positive ^ 1, thus k is 
of hyperbolic type and we have K{ai, 0) — (ai, 0). 

Denote by i>j the eigenvalues > 1 of DK(ai, 0) numerated such that 1 < vi < V2 ^ ■ ■ ■ ^ t^n 

^ /2 

and put 60 = IliLii^j • For a G N we put Ka = boL'~°'. (It is possible that Ka = Ka' for 
a 7^ a'). For any value of Ka and aa,j^i G C, we define 

00 

(9) A; (a, k) = Aa,fc + ^ aaj.Aaf/^"' with al <eE, ^ = 1, . . . , a 

which correspond to the asymptotic expansion for the eigenvalues of an a x a matrix, where 

(10) a = Card {a' I K^' = K^}- 
Let pi be the multiplicity of A; (a, k) as asymptotic eigenvalue. 

Theorem 1 (Gerard, n + 1 odd). For all A > 0, there exists C > such that, if Xa.k is given by 
then there exist coefficients Oa.j^i (z C in asymptotic expansions (0) such that for all N € N 
there exist fcjv G N and cjv G M such that the poles situated in the region 

{Ae C; ImA < A, |A| > C} 

are all in the balls 



A - I Xa,k + / , aa,j,lXaj, 



rriN 



<Cjv|A„,fer^, k>kNeN, 



where is the largest j such that ^ < a is defined in and in each ball there are 

exactly pi poles with multiplicities. 

In his article [S] Gerard reduces the problem to the problem on the boundary of the one 
obstacle, dVli, by introducing the quantum billiard operator M(A) which quantizes the non- 
linear Poincare transform along the trapped trajectory k. The operator M(A) is defined in the 
following way. 
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(11) 



Let i?»,+(A) : C°^ {dn,) ^ C°°(M"+i \ a,), i = 1,2, be the outgoing resolvent of the problem 

(A + \^)H,^+{X)v = m M"+i \ fl. 



extended as an operator i-^ H^^^iCfli) and analytical for ImA < 0. It is known that 

i?i.+ (A) extends analytically as a bounded operator H^^^^dfli) i-^ Hi^^{Zfli) to a domain of the 
form ^ if the boundary dO is C°° and ([3|) if the boundary is real analytic. 
We define Hi{\)v = Hi^+{\)v\QQ^.^^, where d^j, = dVLi and 

(12) Af(A) = il2(A)iIi(A) =7ii?2,+72i?i,+ : H\dni) ^ H\dni), 

where 7^ : -fficlCf^*) ^ H^^^{dn,) is the operator of restriction to Sri^. 

The operator M(A) defined on iJ^((9rii) for ImA < continues analytically as a bounded 
operator H^/'^{dni) ^ H^/'^{dVLi) in the domain of the form (P) if dO is C°° and © if dO is 
real analytic, and there satisfies the following estimate 

3D > 0, 3C> 0, \\M{X)U(HU.^dn,)) < C\\\\'''^'^\ 

where ImA+ — max(ImA,0). 

We have the following relation between the outgoing resolvent i?+(A) in the exterior of O = 
f2i U (u = H+v satisfies ([TT|) with dO, instead of dili) and the billiard operator M, used by 
Gerard in [9], page 91: 
if {vi,V2) e C°°(af^i) X C°°(9172) then 

H+{\){vi,V2) = {Hi^+-H2,+Hi)iI-M)-\i + {H2,+-Hi^+{I-M)-'H2+H2,+Hi{I-M)-^H2)v2. 
Using this relation Gerard proved in [9J that 

Lemma 1. The scattering poles counted with their multiplicities coincide with 

(13) {A, e cr(/ - M{X))}. 
Strictly speaking, Gerard only considers the case of A in a strip 

(14) < ImA < ci, ReA > C2, 

for C2 sufficiently large and n + 1 odd, but his proof also works in more general domains (fTj). ([3|). 

Microlocally near ai, Gerard reduced the problem to the problem of finding the points A 
in a strip ^4]) for which the operator / — e~^"^'^Mo(A) considered in some appropriate space has 
non-trivial kernel. Here A/o(A) is a semi-classical Fourier integral operator associated to the real 
hyperbolic canonical transformation k with a fixed point at (0, 0). Gerard obtained his result by 
approximating of Mq{X) by its linearization. 

In the case when Vj, j = I, . . . ,n are linearly independent over the field of integer numbers 
and the boundaries are analytic, we will use the full Birkhoff normal form of Mq (as in 13]) in 
order to get the explicit formulas for the resonances in any logarithmic neighborhood of the real 
axis of the form ([T]) . 

2. Our result. 

We suppose in addition to ^ that the boundaries Sfii, i = 1,2, are real analytic. Suppose 
that the eigenvalues of the Poincare map 

< t/^^ < . . . < I'V^ < I < 1^1 < ■ ■ ■ < J^n 
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verify the non-resonance condition: 

n 

(15) kj In Vj = 0, kj G Z ki = . . . = kn — 0. 

1 

In order to have good separation of the strings of resonances we impose the Diophantine condition: 

(16) a^(3, \a\,\P\<m | In - /3)| > ^e"^™, > 0. 

Then microlocahy near ai = we can transform A/o(A) into its semiclassical Birkhoff normal 
form to any order r 

Moo(A) :- e-^^F^iIu■■■J.■^/^)^ ^ ^^^q^^ + 1/2) /a. 
Here for any r € Z+, and i in a neighborhood of i = in R", 

(17) 1/A) = ^^o(*) + jF^it) + ±F2{i) + ... + ^F,.(z) 
with Fj polynomial of degree (r — j) in with 

n r 

Fo(i) = G(«,^) + iJ(«), G{i,n)='^^iiti,^j.i=\niyi, H{i) = hj{i) = Oji"^), 

i=l 3=2 

where hj{i) = > d°'H{Q)i" is a homogeneous polynomial of degree j. 
^ — ' a'. 



Here Fq{i) is real when i e I 
We have for any r, r' > j 



[(^) - i^/(l) = O (^jininf '-.'-'l-J + l^ . 



In other words, hj{i) do not depend on r. 

Thus F^,Fj are holomorphic in a neighborhood of i = and Mqo{X) is analytic in any 
logarithmic neighborhood of the real axis of the form 

(18) Aa^b ;= {A e C or A; ImA < AlnReA, ReA > B}. 

Let 

P(A) :=/-e-2*'^^Moo(A). 

As Ijx'^^ — ——{2aj + l)x"^ , the monomials a;" = a;"^ • • . . . • a;"" are formally the eigen- 
2iA 

functions of the operator P(A). We have then 

P(A)a;" = (1 - e-2*'^^i^„(A)) x", if„(A) := e-'^^''(^^i/A) . 
We approximate the problem (|13p with the following model problem: 



(19) 1 - e-^'^^^KaiX) =0 ^ 2dA + Ai^'' (^ ^^^^S 1/A^ 27rfc, a G Z!^, /c e Z, A G A^.b- 

Note that by substituting F"" by its leading term near the origin, ^"^i^T^i^i ■ i-i, we recover the 
pseudopoles ([7]). 
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In general in order to solve we observe that ii F{t;l/\) = Fo (i) + Fi{i) / \ + F2{t) / + . . . 
is either an asymptotic or a finite sum, and F, Fj are holomorphic in a fixed neighborhood of 
1 = 0, and Fo{0) is small (in our case Fo{0) = 0), then 

'2a + 1 



1/A 




dF{'-§^;l/X)2a + l (2gil; 1/A) ' 

di 2i dh 



2d 



for |a|/|A| small. 

We omit r and write equation (|19p in the form 

(20) 

with A/fc as unknown variable. 

It is then clear that equation ([20|l has the solution 




TT 



y =.9(f;^), 5(^;^) = ^+5o(^)+ffi(^)^ + ..., 5o(0) = 0. 
We get the solution in the form 

An TT fli 02 TT -Fi(O) Re (02) / 1 ,„ Im(a2) 

k d k d 2dk k^ \ Adk '^•^ ' fc2 

for k » 1, = 0(|ckp), a — ©(Infc), where aj for j = 1,2,... are polynomials in a (see 
Theorem [5] for details) and -F'i(O) is real constant as the operator M is microlocally unitary (see 



Appendix IA.2[) . The coefficient aj is only dependent on the coefficients in the Birkhoff normal 
form of order j, F^ [i] 1/A) = F^(i) + jF({i) + -^Fl{i) + . . . + ■^F}{i), and independent of r for 
j < r. If Aq and Aq are two such solutions with r < r' then 



k k fc'^+i V 

The main result of the paper is the following theorem: 

Theorem 2. Suppose U5]) and U6]) are satisfied and let Aa^b be the logarithmic neighborhood of 
real axis as in \18\) . For any N Cz N and for all A > there exists B > 0, and r = r(N) € N such 
that there exist functions aj = aj{a), j = 1, . . . , r, polynomial in a of order < j ( aj = O (lap) 
for \a\>l) and there exists bijection b^ between the set of 

and the set of resonances in Aa^b, where elements in each set are counted with their multiplicity, 
such that 

&w(A)-A-0(|Ar^). 

Solutions to (I13|) can be given by means of the Grushin problem for I — M in some weighted 
Sjostrand's space H{AtG) associated to the escape function G. The results of propagation of 
analytic singularities by Lebeau (16' can be applied to localize the analytic singular support 
of M. Outside the analytic singular support the norm of M is small. In a neighborhood of 
the analytic singular support we use the properties of the escape function. Outside some small 
neighborhood Vq of ai € dD,i we get that the operator /— M(A) is invertible. In Vq we reduce M 
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by means of an analytic Fourier integral operator of Bargman type U to the operator e '^^'^^Mq{\) 
in the Birklioff normal form such that for any A^o G 

UM{X)u = e-^''^^MQ{X)Uu + O (|A|-^«) \u\h in H^. 

Here Mq = e-'^^^^'"'/^^ in the sense of |T2]. We use the truncated Birkhoff normal form to the 
order r, which is an analytic Fourier integral operator and the Grushin problem 

in order to show that the solutions of can be approximated by the solutions of the model 
problem 

(21) {A e C, e (7 (/ - e-s^dAg-aF-^c/a/A)^ I 

which leads to ([TO]) . 

Note that as Ro{—X) — RoW it is sufficient to consider the case ReA > 0. 

The structure of the paper is the following: 

In Section [3] wc construct the escape function. 

In Section m we review the method of FBI-Bargman transforms. 

In Section [5] we apply the FBI-Bargman transform to the quantum billiard operator M and 
consider the norm of the transformed operator Mi. 

In Section [6] we reduce Mi to the Birkhoff normal form up to some order r at ai. 

In Section [7] we deform the space near ai such that monomials form an almost orthonormal 
basis in some neighborhoodof ai. 

In Section [8] we derive the asymptotic expansions for the solutions of the model problem (|13p . 

In the last sections CTT^ we prove Theorem [5] by the routine method of Grushin problems. 

In Appendix [X] we review the norm estimates of the quantum billiard operator M in the 
Sobolev space and show that M is microlocally unitary with respect to the flux norm. 

Acknowledgements. The author is grateful to Johannes Sjostrand for numerous discussions 
and constant support during the preparation of the manuscript. 

3. Construction of the global escape function on X = T*drti for the billiard 

CANONICAL transformation K. 

3.1. Definitions of the domains. Let k : T*dfti > T*dili be the billiard map defined in 
dS]). We have K{ai,0) = (ai,0). For p — {x,£,) £ neigh (ai,0) C T*dfli we draw an outgoing 
half-ray 7 issued from p. If 7 intersects 9T*(R"+^ \ ^2) then we denote 7' the reflected half-ray. 
If 7' intersects dT*{M."^-^ \ ^1) then we define /t(p) the projection on T*dili of the point of 
intersection between 7' and dT*{M."+'^ \ Qi). Let D{k) = {p e T*dni; 3k{p) e T*dni} be the 
domain of definition of k. 

We denote X = dfli. The points in T*X can be divided into the following 3 regions: 

H hyperbohc region: {(y,??) G T*X; \rj\ < 1} 

G glancing region: {{y,r]) G T*X; \t]\ — 1} 

S elliptic region: {(y,??) G T*X; \7j\ > 1}. 

We denote pi = (ai,0) G T*X and B{pi,c) := {/? G T*X; dist(/9,pi) < c}. 
We use the following convention: h denotes 1/ReA, A G A^^b- 
Let Wq C T* X he a small neighborhood of pi of the size y^h\n{l/h) given by 

(22) Wo :- B{pi,co^hln{l/h)) C T*dni, 
for some constant Cq > 0. 
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Let Wi C T* X he a larger neighborhood of pi independent of h such that 

Wo CC Wi C ci) C D{k) C H 

for some constant ci > (it wiU be defined later in 

Proposition 1. There exists a real function G £ C°°{T*dQi) with the following properties: 

(1) G^O{h\n\l/h))), \VG\=0[^h\n{l/h)^ , \\/^G\ ^ 0{1); 

(2) for all p G Wo, G{k{p)) - G{p) > G dist{pi, pf; 

(3) for all p€Wi\ Wo, G{k{p)) - G{p) > C/iln(l//i); 

(4) for p e T*dni \ Wi, G{k{p)) - G{p) > Ghln^{l/h). 
Here G is some positive constant. 

We call G a (global) escape function for k. In the following sections we prove Proposition [TJ 

3.2. Construction of the local escape function Gint- As k is of hyperbolic type, by the 
stable manifold theorem, in some neighborhood W of (ai, 0), k has an incoming stable manifold 
A_ and an outgoing stable manifold which are lagrangian manifolds, intersecting transversally 
at (ai,0) characterized by 

A± = {(x,C) e W\ K-"(x,0 G W for aU n G Z±} 

and if (x,^) G A± then k^^^{x,^) i—f (ai,0) exponentially fast as n — > ±oo. 




Figure 1. A part ofT*dili near (ai,0). 

There exists a local canonical transformation / : neigh((ai, 0), T*X) i-^ neigh(0, M^") such 
that /(A+) = = 0}, /(A_) = {x = 0} in a neighborhood W of (ai,0) which corresponds to 
the symplectic change of coordinates. In these new coordinates the differential of k at (0, 0) has 
the form 

D.(0,0)=(^ (AO' 

the eigenvalues of A are z^i , . . . , f„ . 

By the theorem of Lewis and Sternberg (cf. [27], [6]) there exists a Hamiltonian p{x,S,) 
such that K = expHp in the sense of formal Taylor series (cf. [T^] [13]) and p has the form 
p{x,^) = B{x,^)x ■ ^, where 5(0,0) = InA with eigenvalues pj = Ini/j, j = 1, . . . ,n. 

When K is of the form expHp near (0,0) and G G G°° with dG{0,0) = 0, then locally in a 
neighborhood of (0, 0) we get by Taylor expansion 

G(expi/,(p)) = i^P^(P^ ^ ^^P) + HpG + {{x,0') 
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and G{k{p)) — G{p) = HpG + O ((a;,^)^) . Thus G with properties as in Proposition[l]is an escape 
function for p in the usual sense as defined in [lOj. This can be used in order to construct G 
focahy near 0. 

Suppose that the coordinates are chosen such that 

n 

K = expHp, p = ^p,iXi£,i + O {{x,£,)^) . 

i=l 

It is enough to construct an escape function for p := x^, x, ^ e JR. 
Inside the baU + < 0{l)h\n (l/h) we can take 

G^{x,s,)^\{x^^e)- 



Then, HpGo = x"^ + = O {h\ii {\ / h)) , VGo = O (^^/iln(l//i)j , Go = O (/i In (l//i)) . 

For cih\n{l/h) < + < C2 we use the fohowing ansatz: Gi — fh{x) — fh{£.)- 
We have Hp — xdx — £,d^ and we want that HpGi ^ h\n{l/h) in this region. 
We can for example choose fh{x) satisfying 



x dxfh{x) = h\n{l/h), y/h\n{l/h) <\x\<c 

M^h\n{l/h))^h\nil/h)/2 

which gives fh{x) — /iln(l//i) In | — - | H ^ ^ \ And similarly for fhiO- 

\^yh\ii{l/h) J 2 

With the above example in mind we define a local escape function in the following way. 

Let / be a local real canonical transformation f{y,r]) = (a:,^) such that A_ is transformed 

into {a; — 0} and A+ = 0} and let Wi be the image of Wi, 

f{Wi) = t?i C neigh(0,M2"). 

Put s = h\n{l/h) . 

The local escape function Gi on Wi is defined by Gi = G+ — G_ , where 



Then 



G+ = isln(s + a;2) , G_ = In (s + C^) 



„ 1 . /s + x2 
Gi — -sin 



2 \s + e 



li x^ « s then 



G+ = Is Lis) + In ( 1 + ^ U = is In (s) + is ■ ^ + sO i ^ 



2 

We have then for a;^ + < < s 



2x 2^ 



(23) ft = - G_ ^ 1„1„ (1±|1) = - f^) + O + ^ >^ 

If s << a;^ < G, ^ = 0, we get 

Gi = r^'^ '^'^^t^^ ) = r^'-'i^) +^(i)= ^Wl/M)^ + 0(feln(l/Mlnln(l/fe)). 



Thus the variation of Gi is of order /iln^(l//i) 
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For the derivatives we get 



s + 

2x2 



. , ^ s+e' 

As 

^^'^^(^'^) = 7T^ 

the maximum of d^Gi is attained when x'^ = s with the maximum value y/s/2. Using the similar 
estimate for d^Gi, we get 

iVGil = 

We have 

dlG,{x,0^-r^^U-^ 
[s + x'^)^ \ s + a 

and we get that d^Gi attains the minimum at x'^ — 3s with the minimum value —1/8. 
By the similar estimate for 9|Gi, we get 

|V2Gi| = 0(1). 
Similarly we get |a;9^Gi| = 0(s), \£,d^Gi\ = 0{s). 

We denote the billiard canonical transformation in the new coordinates by k := f o k o : 
Wi 1-^ W\. Denote p — (x,^). We have proved the first statement in 

Lemma 2. There exists a real function Gi G C°^{Wi) with the following properties: 

(1) Gi=0 {hln'il/h)) , iVGil = O (v//jln(l//i)) , |V2Gi| = 0{l); 

(2) for all p € Wi, \p\ < c^ h\n(l/h) , Gi{k{p)) - Gi{p) > C |p|2; 

(3) for all p e Wi, \p\ > cy/h\n{l/h), Gi{k{p)) - Gi{p) > Gh\n{l/h). 
Here c, G denote different positive constants. 

Proof of (2) and (3): The canonical transformation in the new coordinates is of the form 
k{x, £,) = {Ax, (A^)^^^) + 0{{x,£,)'^), where A is diagonal with entries e'^\ i — 1, . . . ,n. 
We put n — 1, extension to the general case is straightforward. 
Let s := h\n{l/h). Then, choosing Wi small enough, we have for p e Wi, 

1 /s + t2\ 1 / s + p2^t2 \ 



If |p|2 = .x2 + ^2 < f^g^ then we can use the asymptotics ([23]): 

Gi(k(p)) - Gi(p) = (e2^ - l)x2 - (6-2^" _ 1)^2 ^ (^_Ppj > 2^|^|2 ^ ()Ppj 
If cs < |p|2 = a;2 _|_ ^2 <- jTQj. g gjj^aii enough, we have 
(24) - GM = i. (in + 1.1 (j^^)) + O <»l»'l) S 

using that the right hand side of (j24p is increasing function in x and in ^. Here G^ > depends 
on the size of order e of Wi . | 
In the original space (in coordinates {y,ri)) in a neighborhood Wi C T*dQi we put 

Gint(2/,r/) = Gi(a;,0, /(y, r/) = (x, ^ ^ Gint Gi o /-i. 
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We call Gint the local escape function for k. 

3.3. Construction of the exterior escape function Goxt- As before D{k) denotes the domain 
of definition of k and 

R{k) := {k{p), p e D{p)} 

the image of k. On R{k) we can define the inverse k^^. 
We define 

f2+(0) = H = Hug, h+{l) = Din), for j > 1, := Din^), 

oo 

h+iO) D f2+(l) D f2+(2) D . . . , fl h+U) = A-*, 

i=o 

where ext stands for "extended" by the iterated action of k. Note that as the obstacles are strictly 
convex and p ^ A'^*' then there is m S N such that K™'{p) — ui, infinite point. 

After j reflections lie all reflected rays except may be the last one in some neighborhood of 
(fli, 0). After e-small perturbation of the outgoing ray it will have the same number of reflections 
as an unperturbed one. Thus we have the following lemma. 

Lemma 3. For each j = 0,1, . . . there exists e(j) > such that 

dist{n+ij + i),n\n+ij))>eij)>o. 

Define := \ + 1) C T*X such that the half-rays issued from the points of 

come back exactly j times and then disappears to infinity. For example corresponds 
to the region in Ti. such that if p e ^^+(1) than k(p) — p' ^ T* X but k{p') = ui - infinite point 
(only one reflection). We have 

■.= Hug\D{K), := K~\n+{o)), ... 

OQ 

Lemma[3]shows that approaching the glancing surface (possibly tangentially) f7+(j), j = 0, 1, . . . , 
are well separated. 
We define also 

+ i) {p e n+ij); dist(p, + 1) < dist(p, - 1))}, 
~ i) := {p e dist(p, n+ij + 1)) > dist(p, n+u - 1 ))}. 

Thanks to Lemma [3] these domains extend naturally to £. We denote such extended domains 
again by 

Let be a step function which is constant on each il+(j), j — 0, i, 1, |, . . . , iV, defined as 
follows 

G'+|Uj>ivO+(i) = 0' G'lin+iN-i/i) = -^hlii^il/h). 
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Let uj e C^{W) be such that uj{x) > 0, lu{x) = 0, for |a;| > 1, J^^uj{x)dx = 1. Put 
0Je{x) ~ €^"^00 (x/e) . For u E C"^(M") we consider convolution 

{uJe * u){x) = J uje{y)u{x - y)dy = j uj{y)u{x - ey)dy. 

The mollifier operator u i— s- * u extends to the continuous Hnear operator D' C°° with 

supp(a;e * m) C suppu + {x G M" : |a;| < Ce} 

for some C > 0. 

Using the partition of identity on T* X and moUifiers on any coordinate neighborhood wc can 
construct the regularization G+ G C°°{T*X) such that for p G UjLi ^+0) ^^^d ^ small enough 

G+{K{p))-G+{p)>-h{\n\l/h))>0. 

In the similar way we define G_ associated to ^-{j) defined as before but for the incoming 
half-rays (by replacing n with n^^): 

h^{0) ^H = Hug, = D{k-^) = R{k), forj>l, f2_(j) := D{k-^) = R{k^), 

OQ 

f2_(0) D f2_(l) D f]_(2) D . . . , fl ^-(J) = ^+'' 

i=o 

= {{x,0 e T*X, Vj G N, 3K~={x,i) G T*X}, 

where A™' denotes extension of A+ by the iterated action of k^^. If p ^ A!p* then there is m G N 
such that K~™{p) = UJ. We have an analogue of Lemma [31 
Define := f2_(j) \ + 1). We have 

17_(0) :=Hug\L'(K"^), r2_(l) :=K(r2_(0)), ... 

OO 

+ i) := {p G f^-(j); dist(p,l]_a + 1) < dist(p,fl_(j - 1)}, 

^-U - I) ■■= {P e ^-U); clist(p,0_(j + 1)) > dist(p,r!_(j - 1))}. 

We denote the extension of ^^{j) to £ again by 
Let be the step function: 

G'-|u,>„n+-(j) = 0' G'-|o_(JV-i/2) = 
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We define the regularization G_ G C°°{T*X) such that 



Af-l 

^ G_(a^(p')) - G_(p') < --/iln2(l//j) < 0, p' = K-\p) e U l^-O')- 



Let Gcxt G+ — G_. 
We define 

N N-1 

(25) W2 = W2{N) := y U |J l)-(i). 

i=o 

Lemma 4. i^or aZZ p G W2(A^) we /laue 

Gext(K(p)) - Ge.t(p) > i/lhl2(l//i). 

Proof: If p e i G [1, N], then 3 j G [0, iV - 1] such that p e ^-{j), then we have 

G^^ip)) ~ Ge.t(p) - (G+(^(p)) - G+(p)) - (G_(a^(p)) - G_(p)) > 
>-h\^\\lh)-{--h\^\\lh)). 

I 

Thus we gain at least ift,ln^(l//i) on W2. 

We have that Ai'^' U U^Ar+i^+(*) = ^^+(^ + 1) is a neighborhood of A^.'^S and A^''* U 
\}%N ^-(j) = ^^{N) is a neighborhood of A^'^*. 
Let 

(26) Wi=n+{N + l)r\Q-{N). 

If 

00 

pGAl'^^U IJ 17+(i)\ W^i = f^+(iV + l)\ (f^+(A^ + l)nf}_(7V)) =neigh(Ai'^*)\W^i 

i=N+l 

then G+{p) — and there is j such that p £ Uj!jo^ ^-0) and Gext(K(p))— Gext(p) = — G_(k(p)) + 
G_(p) > \hln^{l/h). Similar for p G A^''* U U°1a, ^^-(j) \ Wi. Thus we gain there too at least 
lh\n\l/h). 

Together with Lemma |3]we get that on H \ Wi we have Gcxt('«(p)) — Gcxt(p) > ■|/iln^(l//i). 
Inside Wi we use Gint Gi o/^^, the interior escape function defined in the previous section, 
and apply Lemma [2l 

Let G be a large constant and x ^ Cg^{T*X) such that x = 1 in Wi. We define the global 
escape function 

G(p) - ^x{p)GnAp) + Gext(p) e C^{T*X). 

Then 

G('^(P)) - ^xi<p))GMp)) + g^Mp)) 

and G{k{p)) — G{p) verifies the inequalities stated in Lemma[TJ | 
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4. The FBI transform of Helffer and Sjostrand and the Bargman transform 

WITH the global CHOICE OF PHASE. 

In this section we remind the definition of the space H{AtG) as in [10]. We follow the presen- 
tation in [23] . 

4.1. Complex domains and symbol classes: generalities. Let X denote a complex neigh- 
borhood of X. Let A e T*X be a closed I-Lagrangian manifold which is close to T*X in the 
C°°-sense and which coincides with this set outside a compact set. Recall that "I-Lagrangian" 
means Lagrangian for the real symplectic form Im a, where cr — da^^ A da^j is the stan- 
dard complex symplectic form. This means that if we choose (analytic) coordinates y in X and 
let {y,ri) be the corresponding canonical coordinates on T*X and T*X, then A is of the form 
{{VtII) + *-f^tG(y,r;); {i/tV) ^ T* X} for somc real-valued smooth function tG{y,ri) which is close 
to in the C°°-sense and has compact support in rj. 

Here Hq denotes the Hamilton field of G. Since A is close to T*X, it is also R-symplectic in 
the sense that the restriction to A of Re cr is non-degenerate. We say that A is an IR-manifold. 

A smooth function a{x,£^; h), defined on A or on a suitable neighborhood of T*X in T*X is 
said to be of class 5'™''=, if 9P5|a = O(l)/i-'"(0'="l«l. 

A formal classical symbol a G S^''' is of the form a ~ h~™{ao + hai -I- . . .) where Uj €E 5'°'*^"-' 
is independent of h. Here and in the following, we let < /i < ho for some sufficiently small 
ho > 0. When the domain of definition is real or equal to A, we can find a realization of a in 
g7n,k (denoted by the same letter a) so that 

N 


When the domain of definition is a complex domain, we say that a G S^''' is a formal classical 
analytic symbol (a e S^J^^) if aj are holomorphic and satisfy \aj\ < CoC^ . 

It is then standard, that we can find a realization a G 5"".'= (denoted by the same letter a) 
such that 



d^d'^d,,^a = 0{l) 



J2 h^a^ 



0<J<l(C)l/Co/i 



< 0(l)e-l<«>l/c^l^ 



where in the last estimate Co > is sufficiently large and C, Ci > depend on Cq. We will denote 
by S^'^ and also the classes of realizations of classical symbols. We say that a classical 

(analytic) symbol a ~ h~"^{ao + hai + . . .) is elliptic, if ao is elliptic, so that £ 5°'"*^. 

4.2. The FBI transforms. Let X ~ dfli. Following Sjostrand and Zworski [53] we intro- 
duce an FBI transform which is a map T"^ : C°^{X) ^ C°^{T*X) given by 

(27) T"^u{a,h)^ f e"t'^°'''-'^'^a{a,x]h)x{a,x)u{x)dx, 

Jx 

where x is smooth cut-off with support close to the diagonal 

A = {{a, x) e T*X X X, a^= x} 

and equal to 1 in a neighborhood of A, a g S^'^ is elliptic, and cj) is an admissible phase 
function. There exists b{a,x;h) G S^'^ such that if 

(28) S"^v{x,K)^ ( e-"'''^°'''='>/^b{a,x;h)x{a,x)v{x)dx, (l)*{a,x) = (j){x,a), 

Jt'X 
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then S^^T^^u = u + Ru, where R has a distributional kernel R{x,y; h) satisfying 

Let G e C°°{T*X) be the escape function constructed in the previous sections. We define its 
C°° extension to T*X, the complex neighborhood of T*X, with the property that 

{dG)p\jT^^T'X) = 0, 

where J : Tp{T*X) t-^ Tp{T*X) is the complex involution. Denoting the extension by the same 
symbol, G G C°° {T* X), we now define a C°° /— lagrangian, i?— symplectic submanifold of T*X : 

(29) AtG = exp{tH}f"'){T*X) C T^, lma{x, Hji'"') = dG{x), cr = da^ A da,. 

Here by a C°° manifold we mean a graph of a C°° function. We remark that a C°° graph 
in T*X over T*X which is /—lagrangian (i.e. its almost everywhere defined tangent plane is 
Lagrangian with respect to Ima) can be locally written as a graph of a differential of a C°° 
function, G onT*X ~ T*W' (locally): 

dG dG 

{a = (aa;. Off) e r*C" : Ima^, = — {Reax,Reac), Imae = (Re a^;, Re at )}, 

oRea^ aRea;^ 

(Rea:r,Rea^) G r*M" ~ T*X, G = tG. 

The form —lm(^dz\\^^ is (formally) closed and, as Ajg is close to T*X, it has a primitive which 
is a G°° function on Ata ■ We denote it by i/ G G°° (Ate ; R) . We normalize H by demanding that 
it is equal to for |a^| large enough. 
Parametrizing Ate hy T*X : 

AtG^{ix + it—,C-it—); (x,e)Gr*X}, 

we get 

(30) Htia,,a^) - • td^G{x, + tG{x, 0, («.,««) = i^, + ^tHaix, 0, (a^, e T*X. 

The weight function Ht{ax, a^) have the same properties as tG due to Proposition [TJ Thanks 
to the global properties of the amplitude and the phase function of T^^ we can continue it in a 
to a neighborhood of T*X, T*X, in T*X and in particular we can define 

A deformation argument as in Sjostrand gives an approximate inverse S\^^ , with the same 

properties as S above, defined by ([28)1 with T*X replaced by Ate- We put for u G C°° 

(31) 

ll"ll?G - Ikll W,(«,>") = rA*o^^lli.(A.o.(«,>'"e-*/'-) = / |TA,,^(«;Mn(a^)p-e-2^'(")/''rfa. 

■J AtG 

Using S'ajq we can show the independence of this norm of the choice of a specific phase function. 
From 23J and 29J it follows also that the kernel of Ta^q SAtc ^^le form 

TmoSa.o («, P) = c(a, /3)e-t'A("./') + r{a, f3; h), 

where ip{a, 13) = v.c.z{(t>{a, z) — (f>* ((3, z)), c(a, P;h) = 0{h^") is supported near a = /3, 

r(a,/?;/i)=OAr(l)e-VC/«|(/3c)r^ 

and 

-H{a) ~ lmil;{a, P) + H{/3) < -^-d(a,f3)^. 

c 
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Here d{a,(3) is any non-degenerate distance on Ate- We used that |V^G| = 0(1). 

4.3. Generalities on the Bargman transforms with the global choice of phase and 
choice of the norm. Following [23] we replace the norm || • ||fG by an equivalent norm obtained 
by decomposition of T*X. 

In the region where ^ is bounded, it will be convenient to work with transforms which are 
holomorphic up to exponentially small errors, and for that we make a different choice of T, and 
take an FBI-transform of Bargman type with a global choice of phase. 

The standard Bargman transform in M" is given by 

Tau{z,h) = h-^ I e-^'-''^^/^*'u{x)dx, Tq : ^^(M") ^ H^oiC), $o(^) = (Im2)V2, 



where iJ$ := {u holomorphic, J \u\'^e ^"^^^dx < oo}, associated to the complex canonical trans- 
formation 



Let (f>o{a, x) = i(Im z)^/2+i(z— x)^/2, z — a^^ia^. Then (l)o{a, x) is admissible phase function 
and thus e~^*"*^^^/'*To is local FBI transform of the type described in the previous section: 



KTo(AtG) = A*,, = $o(z) +<G(Rez,-Imz) + C'(t2). 

We need the global version of Bargman transform on T*X. We equip X with some analytic 
Riemannian metric so that we have a distance dist and a volume density dy. 

The function dist(x,?/)^ is analytic in a neighborhood of the diagonal in X x X, so we can 
consider it as a holomorphic function in a region 

{x,y) e X X X; dist(x,y) < ^, |Ima;|, |Ini?;| < ^ 

Put 

(32) (j){x, y) = iAodist(x, yY , 

where Aq > is a constant that we choose large enough, depending on the size of the neighborhood 
of the zero section in T*X, that we wish to cover. 
For X e X, |Imx| < 1/C, put 

(33) Tu{x;h)^h^ [ ei*^'''y'>x{x,y)u{y)dy, ueV'{X), 



where x is a smooth cut-off function with support in 

{{x,y) eXxX; \lmx\ < 1/C, dist(y,y(x)) < 1/C} 

containing HU G- Here y{x) g X is the point close to x, where X 3 y i-^ —lm(f){x, y) attains its 
non-degenerate maximum. We collect the following facts from [23j : 

• The function ^a{x) = —lm(fi{x,y{x)), x X, |Imx| < 1/C, is strictly plurisubharmonic 
and is of the order of magnitude jlmxp. 

• A$Q := {{x,id^o) G T*X} is an IR-manifold given by A^^ = kt{T*X), where kt is 
the complex canonical transform associated to T, given by (y, —(j)'y{x, y)) !—> (x, 0^(x, y)). 

Here and in the following, we identify X with its intersection with a tubular neighborhood 
of X which is independent of the choice of Aq in 
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• If = L^{X; e-'^'^°/^L{dx)), for L{dx) denoting a choice of Lebes gue measure (up to 
a non-vanishing continuous factor), then T = 0(1) : L'^{X) ^ L|^, 'd^T = ©(e^^/^'') : 
L'^{X) I— L|, . This means that up to an exponentially small error Tu is holomorphic 
for u e L^iX) (and even for u e V'{X)). 

• Unitarity: Modulo exponentially small errors and microlocally, T is unitary 

L^{X) ^ L^{X;aoe-^'^°^''L{dx)), 

where L{dx) is chosen as indicated above, and ao(a;;ft-) is a positive elliptic analytic 
symbol of order 0. 

• Let htG C T*X be an IR-manifold as before. Then KT^h^ta) = A$t7 where $t can be 
normalized by the requirement that <i>t = $0 near the boundary of X. (Here is where 
we have to choose Aq large enough, depending on Ktc- In the applications, for a given 
elliptic operator, Ktc and T* X will coincide outside a fixed compact neighborhood of 
the zero section, and the whole study will be carried out with a fixed Aq.) 

Locally we can identify T*X ~ T*IR" and take (j){x,y) — i{x — y)^/2. Then locally kt : 
(2/: V) ^ [y ~ v) =• '^ith <i>o — (Ima;)^/2. Using the local formula 

-ff* (2/; ^7) = -R-e 77 • Ini y + tG{Re y, Re r?) 

we get 

^t{x) = 'y-c.(y,Re^)eC"xR" (-Im0(a:,?/) - RC77 • Imy + tG{Rey,RcT])) ^ 
Mx)^Mx)+tG{x)+Oif\xf), 

(34) where Gix) = Giyix),r,ix)), {y{x), r^{x)) = k^' (^x, 

Let be a neighborhood of pi e T*X, such that Wi CC W. We identify W with nriW). 
It follows from [10], Section 9, that for u £ H{A; (a^)™) we have 

\\Tu\\h^^{W) < c\\u\\tG- 

Let W CC W and Wi CC W be another neighborhood of ai which is identified with {a G 
AtG', Rea G VF'}. One can show, using methods from [TU], [12], that for < i < io with to small 
enough i3 > such that 

\\u\\tGM' < {\\Tu\\H,^iw) + e-'''/'^\\u\\,a,Cw' 

where in general the seminorm ||u||tG^vK is defined as in (|3ip by integration over W C Ate only. 
We define a new norm uniformly equivalent to || • ||tG as ft- ^ by 

hlltG.CM/' + \\Tu\\h^,^^w)- 

5. The FBI-Bargman transform of the quantum billiard operator M. 

The operator M was defined in the Introduction, formula p2p . Here we give an alternative 
definition (see Burq [4]). 

Let h = 1/A G M4.. If Hi,^ is the outgoing resolvent of the problem 

{-h^A~ l)Hi^+u =0, inCl^i 

we consider the outgoing resolvent of the problem 

{-h^A - l)i?2iu = 0, in CrJa 



(35) 
we c( 
(36) 
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The operator M is defined as 

(37) Mu = H2iu\QQ,^. 

Some known facts about Al are collected in Appendix [Al 

Let Wi be the neighborhood of pi = (ai,0) G T*X independent of h and defined in (|26p : 
Wi = Q,j^{N + 1) n ri_(A^). Let G be the global escape function as in Proposition [TJ 

5.1. Properties of the kernel of M on ZWi. Let T be the FBI transform of Helffer and 
Sjostrand as in (j27p with admissible phase function (/> in a complex neighborhood ?7 x of {pi, ai) 
and let S be an approximate inverse to T as in (|28p . Let fJ-, p € T*X be arbitrary real for the 
moment. As M is bounded operator L^(X) i-^ i^(^) and T, 5 are bounded L'^iX) i-> L'^{T*X) 
respectively L?{T*X) i^(^) then the distributional kernel K of TMS* given by 

(38) TMSTu{a,h)= f K{a,f3;h)Tu{f3)d[3, a e neigh{pi), 

J neigh (pi) 

is well defined. 

By the results on the propagation of Gevrey 3 singularities due to Lebeau [16] (see also Burq 
[1]) it is known that on the complement of the graph (k) : 

(39) e T*X X T*X, a = k(/3)}, 
K{a, f3; h) is small: 

(40) K(a,^;/i) =C'Ar(l)e-i/"'*'''(max(|(a^)|,|(/3^)|))"^ iia^K{l3), for any iV > 0, 

if k{p) is defined, otherwise no condition. Since H{a) has compact support and is small, we have 
the same estimate for the reduced kernel: e'^^^'^"^+^'^'^^^/''A'(a, /3; h). 

In a neighborhood U xU oi the graph of k ([39]) . from the general principles, as M is bounded, 
we know only that K{a,(3;h) is of order 0{h^'^ ) for some N' > 0. If a,/? G CWi we can get 
a better bound on the reduced kernel e(^^(")+^(^))/''A'(a, h) by using the properties of the 
escape function G, LemmaHJ if a = k(/3) e T*X \ Wi then G{a) - G{(3) > Gh\n{l/h) and the 
bound is still valid in the whole neighborhood of {{a, f3) G T*X \ Wi; a = k{(3)}. 

Let AtG be the IR-submanifold of the complex neighborhood of T*X defined in ([^5]) . Let 
K{a, f3; h) denote also the extension of ([55)) to the complex neighborhood of {pi, pi) of the form 
{U n Ate) X ([/ n Ate)- Let H € C°°(AtG;R) satisfy daH^/^^^ — —lma^dax\KtG ^'^^ t)e equal to 
zero for \a^\ large enough. By ([SO]) we have locally 

Ht{ax, ct^) — — Rea^ • Ima^ + tG(Re a^:, Rea^). 

Then for {a, (3) in a neighborhood of {(a, f3) G CW^i; Re (a) — K{Re /?)} we can arrange that 

(41) e-^(")/''K(a,/?;/i)e^(«/'' =e-ckmin(/.h.^(iA),|a-K(;3)|^) (naax(|(a^)|, K/S^)!))"^ 

for some A^o > and for any A^ > 0. 

Then by the general theory of [TUj and [53] the estimate ([¥T]l implies that for u G H{AtG) we 
have the bound: 

(42) \\Mu\Ua,Cw' < 0{h^")\\u\UG 

for some A'o > 0. Here W is any /i— independent neighborhood of pi, W\ CC W CC W 
identified with {a G Ate! Rea S VF'} and in general the seminorm IjwHtG.H' is defined as in ()3ip 
by integration over W C Ate only (see the end of Section [473| . Boimd ([42]) extends to complex 
^ S B with some different power > 0. 
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5.2. Local form of M. Microlocally near (ai,0) G T*dVli, the asymptotic solution to ([55]) is 
given by means of WKB construction. 

In the hyperbolic zone 7i, M is approximated by Fourier integral operator H e l'^{X, X] k!), 
where k' — {{x,^;y,—r]) : (x,^) = K{y,r])}, associated to the real canonical transformation of 
billiard k, with the real phase and can be taken in the form (see [9]): 

(43) Hu{u, A) = (^^^ J g-a(.(x,e)-,e+2d)^(^^ Q. x)u{y)dyd6, 

where s solves the cikonal equation, s(a;, 9) — yO parameterizes the canonical transformation of 
bilhard k, s(0, 0) = 0. 

If A is complex, A e ^a.b^ then we denote Ai = Re A, A2 ~ ImA and write 

Hu{x, A) = J I e-^>^^«-,o)^yO+2d)l^^^ q. ^^^^(^^ y^^^y^Q 

with 

6(x, y, 61; A) (1 + z tanarg e^^'^''^^^-f>^-y^+'^%{x, y, 9; Ai). 

Let T : L'^{X) H^„, H{At) H^^ (modulo exponentially small errors) be as in ([33|l 
associated to the complex canonical transformation kt such that kt{T* X) = A$q and kt(A4g) i-^- 

For <& denoting either $0 or <I>t, the Fourier integral operator H : H<j, 1-^ is associated to 
the real canonical transformation k in the sense that locally k(A$) ~ A$. 

Let Vi = tTxKtWi and for any x G Vi let 7(x) € Ate be such that TTxHxi'^ix)) — x. In 
Wi C 7i we apply the Bargman transform and denote H the transformed operator associated to 
the canonical transformation k — kt o n o , satisfying for u G i^(Wi), 

TMu = HTu + 0(/i°°)||m|| in i/*. 

As ffw is asymptotic solution to ([55]) and obstacles are analytic then _ff can be taken analytic 
of the form 

(44) Hu{x,\)^ (^-^^ J j e-'^'f>^'''y^'^^a{x,y,0;X)u{y,X)d9dy, u€H^. 

with analytic phase and an amplitude which is an analytic symbol of order realized with some 
suitable contour T{x) C Wi passing through the critical point, and introducing some cut-off 
functions. 

It is well-known that, by choosing an appropriate integration contour, H is associated with a 
kernel G ^^#(2:)+$(!7) such that, formally: 

Hu{x,h)^ [ Kfj{x,y,h)u{y,h)e~'^'^^y^''''L{dy), u&H^. 



The kernel is uniquely defined as an element of -ff$(i:)+$(27) through the data of H, but not 
as a function. 
Denote 

k{x,y;h)=Kfy{x,y,h)e-^''^yy\ 

If Re 7(0;) 7^ K(Re7(y)) ^ Wq then by the results on the propagation of analytic singularities of 
Lebeau [TB] (see also Burq i4 ) we have 

(45) e-*(")/H(a;, y; /i)e*('')/'' = O (e^'^'') , 
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where $ is either $o or <i>t = ^o{x)+tG{x) + 0{t'^\x\'^), where G{x) := G'(Re7(a;)), nxKrilix)) = 
X. Here we used that G = 0{hln{l/h)). 

If Re 7(2:) = K(Re7(i/)) then we use the properties of the escape function: 

-G'(K(Re7(2/))) + G(Re7(2/)) < -C/iln(l//i), c> 0, if RC7 e VFi \ VKq. 

Choosing Wq large enough (cq in the definition is sufficiently large) we can arrange that 

(46) e^*'(^)/''fc(x,?/;ft)e**(^)/'^ ^^^(/i^") if Re7(a;) = K(Re7(y)) n - 0, 
for some A'o > and Vq = tTxKtWq. 

5.3. The Bargman transform reducing the stable manifolds to a; = 0, ^ = 0. The 

outgoing A+, incoming A_ stable manifolds for k are lagrangian, intersecting transversally at 
(0,0). We can introduce real symplectic coordinates {x,£,) such that A+ : ^ = 0, A_ : x = 0. 

We consider the image of A± under the complex canonical transformation kt '■ T*X 1-^ A$q 
associated to the Bargman transform T. Transformation kt preserves the properties of A±. We 
will often write A± instead of kt(A±). 

We know that the fiber {x = 0} is a Lagrangian manifold, which is strictly negative with 
respect to A$q. It follows that {a; — 0} and k(A„) are transversal. Then A± are given by 
^ — dx4'±, where (j)± are holomorphic, and — Im0+ — ^o{x) < 0. The image of Ate is of the 
form A$j, where $t(x) — $o(a;) +tG{x) +0(^^1x1^), where G is considered also as a function on 
A<[,„ (see (121): 

(47) G{x) := G o {Tr^ o kt|t*x)"^ = G{y, r/) with {y, 77) e Ate given by Tix{>iT{v, v)) = x- 

(If <j) — i{x — 2/)^/2, and {y, rj) e K^" (instead of complex Ata), then x — y — if] and G{x) = 
G(Rex, —Imx).) 

The strict positivity of A+ with respect to Ate then implies that 

Im0+ + ^t^ 

Using the strict negativity of A_ with respect to Ate we have (see [25] ) 
Lemma 5. There is totally real linear space, L C C", of real dimension n, such that, 

$t + lm(f>_ ^ -\x\'^ on L. 

Considering still the situation after application of kt, we now introduce a new canonical 
transformation Kp which maps A+ to — 0} and A_ to {x = 0}. Then it is easy to see that 
Kp is given by kf ■ (y, — 0'y(a;, y)) (x, 0^(x, y)), where the generating function (/)(x, y) verifies: 
det<^^',j,^0, 0(O,y) = -0_(y). 

Let /(x, y; h) be classical analytic symbol defined near (x, y) — (0, 0). Using the Lemma[Sl we 
see that if it S H^^ , then we can define Fu G , by choosing a nice contour for the integral 
expression, 

Fu{x-h)= j e'*^^^y^'^f{x,y-h)u{y)dy. 

Here, we use terminology of |22j . is a new strictly pl.s.h. function, determined up to a constant 
by the relation A^,^ = kf(A$J. 

As A+ is strictly positive with respect to A$j we get that Ki?(A+) is strictly positive with 
respect to A|,^ and the strict positivity of hf{Aj^) : ^ = (that we shall from now on denote 
simply by "A+") means that. 
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Up to exponentially small errors (modulo equivalence in the spaces H^^^) we can invert 
F by an operator, 

Gv{y-h)= [ e-''*'^''^y'''^g{x,y]h)v{x)dx. 



The canonical transformation of billiard k is transformed to 

K = exp Hp, Dk(0, 0) [A^y^ ' ^ " ^ diag(/ii, . . . , ^„), 

where we have used the Lewis-Sternberg theorem (see [13 j) and = denote the equivalence relation 
for formal Taylor series at (0,0). 

We denote Wi := kf(W^z) and n{W.,) for i = 0, 1. 

Composing the Bargman transform T with the integral transform F we get an operator of 
norm C(l), 

(48) FT: H{KtG)^ H^{Vi). 

Here i > is small and fixed, G is the escape function introduced in Section [3l $ denotes a 
function having all the properties of and Vi is small open neighborhood of in C". 

The direct definition of T and F only gives that d{FT) is exponentially small, but we can 
correct this by solving a 9-problem, using the fact that $ is strictly plurisubharmonic. 

Composed operator FT has microfocal inverse SG of norm 0(1) : i?$(Vi) ^ H^Ata) with 
the properties: 

• SGFT is a pseudodifferential operator of order adopted to A^g in the sense of [lO] . 
which has compactly supported symbol and which realizes the identity microlocally near 
pi = (ai,0). 

• We have \\FTSGu — ^Wh/^CVi) ~ ^(^°")ll^llir4.(Vi)7 where Vi CC Vi has the same prop- 
erties as Vi. 

We have 

Theorem 3. The transformed operator FTHSG = e^'^^'^Mi is an analytic Fourier integral 
operator given formally by 



(49) 



Mim(x, A) = j e-'^'^'^^-''^~y'Mx, y, 6; X)u{y, X)dyd9. 



Here b G S°^^, b{0,0;X) = (i^i ■ . . . ■ Vn)"^^^ + 0{\X\-^)., where 1^1 , . . . , f „ are the eigenvalues 
>1 o/Dk(0,0). We have ip{x,9) = A^^x ■ 6 + 0{\{x,9)\^), where A is diagonal with eigenvalues 

Vl,..., I^n, Vi = e^*. 

Ml can be realized as bounded operator H,^{Qi) H,^(fl2) with domains CC CC Vi, 
where we identify fii with ttk (A$ n 7r~^r2i) . 

For any x € C^(C") with suppx C Vl \ ^ and some Nq > 0, we have 

(50) llxMiii|L|(o,)=0(/i'^")||ii||L|(no- 

The form of Mi in (gS]) follows as in The estimate ^SU^ follows from (gS]) and 
Let $t denote the new <tt after the transformation F. 

6. SeMICLASSICAL quantum BiRKHOFF NORMAL FORM. 

In this section we perform an analytic Birkhoff transform k,. up to some fix order r, 

Kr ■■ A$t i-^ A|.^ . 
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Let A € I^A,B and we put h = 1/Re A. We need the following notion of equivalence used in papers 
[13], [Hj. Here p denotes either pair (x,^) or triple {x,y,9). 

Definition 2. Let U E /"(X, X, k') and U E l'^{X,X,k') be two Fourier integral operators. 
Then U = U to the order r if k{0, 0) = k(0, 0) — (0, 0), k, k agree to the order p^''"*'^ at (0, 0), the 
terms with number j, with < j < r, in the asymptotic expansions (in powers of h, corresponding 
to h^ ) of (j),4',a,a (phase, respectively, amplitudes) agree to the order p'^^''-3^+^ at (0,0,0). 

We write a = 0(|pp'"+^ + /i''+^) to denote that a is equivalent to zero to the order r, and use 
= to denote the equivalence to infinite order. 

The theorem proven in 13J for real hyperbolic k says the following: 

Tiieorem 4. Let Mi be a Fourier integral operator as in Theorem\^ which quantizes an analytic 
canonical transformation k : neigh (0,C^") — > neigh (0, C^"). Assume that the eigenvalues of 
dK{0) satisfy 

< I/^^ < . . . < < 1 < 1^1 < . . . < l^n- 

Let pj — log . 

Then there exists a pseudodifferential operator P{x,X^^Dx',h) with symbol P{p;h) ^ pifi) + 
hpi (p) + . . . , such that 

(51) Ml = e~'^^. 

P is uniquely determined modulo '= " and up to an integer multiple of ^ by ^51]) and the choice 
ofpo such thatp{p) - po{p) = 0{\p\^). 
Suppose 

n 

(52) kjpij = 0, fcj e Z =^ /ci = . . . = fc„ = 0. 

1 

Then there exists elliptic Fourier integral operator B and a classical symbol of order 0, 
F{i- 1/A) = Fq{i) + \-^Fi{i) + X-^F2{i) + ..., (h, . . .,in), ^J = ^^2;^, 

71 

Fo{^)^J2^'J'^+^('^^ i?(*)=0(z2), 

such that 

P = B-^F{I;1/X)B, Ml = B-ie-*^^(^'i/^)B, / = (/i, . . . , /„), = (x, a., + 9,^- ) . 

When All is microlocally unitary near (0,0), F{I]l/X) can be chosen to be microlocally self- 
adjoint and unitary, respectively, at (0,0,0). 

Fix some order r > 1 and let 

F'-(z; h) := F(,{i) + X-^Fl-\i) + X-^F^-^{i) + . . . + A-^F^, 
where f J""* is polynomial of degree (r — j) in i {F^ is constant) such that 

F,i^)~Fp{^) = Oi^^'^+'), 

{Fj^-' -x) is resonant up to order 2(r — j)). For any r' > r the terms in expansion of F^ (i; 1/A) 

and F^ with the same indexes coincide, thus the coefficients in expansion F'^(i; 1/A) are inde- 
pendent of r. 
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Then F''(/;l/A) is analytic pseudodifferential operator. From the construction of Birkhoff 
transform B for the Hamiltonian it is easy to see that there exists an analytic Fourier integral 
operator Br associated to an analytic Kr such that 

P-B-^F-{I;l/X)Br^Q'^, 

where is pseudodifferential operator with the symbol q^{p; 1/A) ~ qo{p) + X^^qi{p) + . . . , with 
qj{p) = + /I'-J+i), if < j < r and qj = Pj for j > r + 1. Here B-^ is analytic 

Fourier integral operator satisfying B~^Br — I — and Q"^ of the same type as Qr- 
Operator e~^^^ ii;^/^) jg analytic Fourier integral operator which satisfies 

We denote 

Afo := BrMiB-\ Mo - e''^''^^^'^'^\ 

Operators Mq and i?'' have the form (|49)) . Let 

= A- //e--'<-..«.6',..,..,W« 

and 6 ~ '^''^i' = C'((|a;| + |?;| + |6'|)^'^''^^'++^). The integration contour is chosen such that 

+ lm^'-(x,y,0) + $(2/) ~ -[x^ - \y\^ - \e\\ 

We have 

(53) e^^'^R'e'^'^ = 0{\p\''''+^ + h^+^). 

Relation extends to A e Ka.b with the bound O (h-^ (\p\^''+'^ + h'+'^)) for some C > and 
usual convention h — 1/Ai, Ai = ReA. 

Using (f53|) and that for A S ^a,b, we have le"*^"*^! < C'(1)A™ for some (positive) m, we get 

Lemma 6. Let Wq ~ 5(0, cq-v/s), s = h\ii{l/h), co > as m llg)] and Vb swc/i i/iat 7r(Wo) = Vq. 
Suppose A G Ayi^B. For any G N f/iere is r = r{N) G N smc/i that 

where 0,2 CC fii CC Vq and $ is either $o or $j or any pl.s.h. function close to $o C*^- 

Proof: In the proof h — 1/X is real, extension to A G A^.s is straightforward. 
We write p — {x,y,6), (f){p) — —[tf(x,9) — yO). Stationary point of with respect to y,6 is 
given by (y, 9) — (0, 0). If Mqu = w we have, 

(54) v{x) = /i-" / / ^^{^Ap)+^Ap))/h (^r^^^ ^ ^r(^^ u{y)dyde, 

J Jr(x) 

with a good integration contour r(a;) passing through the critical point and ^r{p) ~ 0{\p\^^^^), 
e ~ with, for < j < r, c] = Cl(|p|2('-J)+i). If \p\ < 0{s^/^) we get 

A{p)=0{s^+'/'), c^{p) = 0{s^+'/'). 

Let C B{xq, C2S^^^), be a neighborhood of xq ~ 0, and let 

ni = 7r3^K"i(7r^if22 n A*) CC B{yQ, cis^/^), yo = 0. 

By modifying Mqu with 0{h^°) in L|(f22), with some Nq > 0, we can restrict the integration 
contour T{x) in ([54]) to |y| + \9\ < es^/^, with e small enough. 

It is clear that e~^^^ (^^i/-*') ig -virell defined as a bounded operator H^{^li) H^{fl2)- 
We want to estimate the difference Mqu — e~^^^' ii'^/^)u on Q2- 
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We make substitution x = s^^"^!:, y — s^^^y, 6 = and write u{y) — u{y), v{x) — v{x), 

p = (x, y, 6) and s^/'^p — {s^/^x, s^^^y, s^^^O). Then, if AIqu — v we have 

iii) ^ (^^^ 'Y I e^i*As'''p)+o{s''^''')) (b^[s^/^p^X) + 0{s--+^'^)) u{y)dyd0 

Let r be the image of T{x) by substitution. We have then x £ C B{0, C2), and for {y, 6) £T 
we have \y\ + \0\ < e. We write (j)'^ := ■i(/)r(s^/2^). Then on f 



Then, if 



-S-^^is^/^i) - Im0^ + S-l$(sl/2£) ^ _ (^|5|2 + |y|2 ^ |^|2^ ^ 

(55) vrii)=(^^^ je^*'^^bl{s^/^i,s^/^y,s^/^6)u{y)dyde, 

we get 

where ||.||$,r2 denotes the norm in H,j,{il). Here > can be chosen arbitrary large by choosing 
r large enough. 



7. Deformation of the escape function in a neighborhood of (0, 0) and final 

definition of the space. 

We write h — 1/Ai, Ai — Re A. Let Wq be the neighborhood of of the size cq as 
in ([22]) (with appropriate identifications of the domains). Let Wi DD Wq be the A— independent 
domain as in ([26]) . 

Operators Mq = B^MiB^i, e^^-^-^'t^'i/^) = Mq - i?'' can be extended to the whole A- 
independent neighborhood of 0. Let U = BrFT, where T is the Bargman transform with a 
global choice of phase ([33| , F is the analytic Fourier integral operator quantizing the symplectic 
change of coordinates near ai = 0, Br - the Birkhoff transform up to order r. On Wi \ Wq we 
can estimate the norm of as in ([50|) . Together with Lemma |6] it implies 

Lemma 7. For any N E N and A e Aa.b there is r E N such that 

UMu = e-*-^2<ig-»AF'-(/;i/A)^^^ ^ e-'^'^'^BJ-u, m L|(Vi), and 
\\e-'^^''R-u\\Li(y^) < Or{h^)\\u\\HiA,a)' where ^{W^) = V^. 

In Section[3]we constructed a global escape function G in r*(9f2i), such that G{k{p)) — G{p) > 
on T*d{^li) and such that we have strict inequality, outside an arbitrary small neighborhood 
of (ai,0) e T*dni. 

Let fiint and be a A-independent neighborhoods of (ai,0) (after usual identification of the 
domains), such that Wq CC CC flint CC Wi. We define a preliminary space Hpro associated 
to the Lagrangian space h-ta as follows 

In some small neighborhood of (0,0), 



(56) Woo = {p, \p\ < cooy/hH^/h)} cc Wo, 

we deform A$j in the following way. 
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We have kjj :— kb^ ° up ^ kt Ate '-^ ^tG- denote again G by G. With G we associate 
an IR-lagrangian manifold which is essentially 

A*G = = e^piitHG)iy,v); {y,v) e T*(neigh(0))} 

for < t < 1, smaU. 

Locally we have A+ = 0}, A- = {x = 0}. For x^ + f < 0{l)hlTi{l/h) we have 

G ^ {x'^ — (^)/2, and Ktc is given by 

Xj — (cost)y.j — i(siTit)7'i.j „ i ^ ■ ^ 

^ ~i{sint)yj + {cost)rij ' '^J' ' - - 

Then KuiAtc) = A*, = U = f $t = i(cotO(Imx)2 + i(tani)(Rex)2. 

Notice that $^/4 = |a;p/2 and that the corresponding IR-manifold is given by ^ = —ix. 
Following [5] we look for a new IR-manifold A which coincides with A^/4 = A$^^^ in some 
neighborhood of (0, 0) for a;^ + < 0{l)h\n{l/h) and with A$j outside. 

First we notice that if qi{x) and q2{x) are strictly convex quadratic forms then we can find 
smooth strictly convex function (j){x), with (l){x) = qi{x) near and with (f>{x) = q2{x) outside 
F, where V is any given neighborhood of 0. Apply this with V = {x 1^ C", \x\ < K} for some 
K = 0{^h\n{l/h)), q,{x) = <f^/4ix), q^ix) = 

Then replace 0(a;) by ^o- = (-^) 0(^), < cr <C 1, in order to decrease the neighborhood of 
0, where 7^ $4, even further while keeping bounded in C^. We have a ^ \Jh\n{\/h) and the 
derivatives of higher order than 2 are diverging as /i ^ 0. Thus we get that the deformation (p^ 
is strictly convex function which is close to the original function in C^. 

The deformation of A^^ is denoted A|,^ and the new weight function satisfies 

= |a;|V2, < e/iln(l//i), ^t{x) = $t, jxp > cqo/i ln(l//i), 

where < e ^ cqo are some constants. 
The final version of the space is given by 

T"^u e Ll(.{T*dni \ W) 



8. Model problem. 

In this section we derive expressions for the coefficients in the expansion in powers of the 
solution of equation ([19]). 

We suppose that fij satisfy the non resonance condition (|52p . Thus all /ij and Vj = e^^ are 
different. We also have fij > ^ I'j > I. 

For A e Aa,b, as in (fT8|) . consider the equation 



We omit r in the notations and proceed in increasing the order of generality. 

A. Assume F{i; 1/A) = G{i, /i) = fM^i-, * — x^, is independent of A and linear in i. 

2dA + AGf-i-(2ai + l),...,-i-(2a„ + l),/i) ^2nk, keZ ^ 
\2iX 2iX J 

2dA + — V Aij(2a, + 1) = 2iTk ^ X(k, a) = fc^ + A V ^J■^{2a^ + 1), (k, a) e Z x W 
2i ^ a 4o 
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We have 

ReA = ^, ImA = ^^j(2aj + 1) => jal^IniA, 



as ReA ^ oo, \a\ = 0(|ImA|), |a| = 0{k^-'-), \a\ = 0(log|A:|), k ~ Re A. 

B. Assume that F{i) = G{i,fi) + H{i), polynomial of degree r, is independent of A, where G 
is as above and H{i) — 0{i^). We have 

r 

(57) =^A-^7.,(y), |1/A|^0, M < lA^-, 

for some e > 0, where \hj{y)\ — 0(|?/p), hj{y) = Yl,\a\=i '^9°'H{0)y" is a homogeneous polyno- 
mial of degree j. Then we have 

A 



where hj{a) — hj (^^7^) = O (lap) is homogeneous polynomial in a. 
It is clear that equation ([58|) has solution in the form 

A TT ai a2 

Then ([55)1 transforms to 

^ + ^ + ... + ^E..(2.. + l) + A^A-,^.(.).0, 

1=1 j=2 

where the last term is of order Odap/fc^). 
Thus we get 



1 

«i = -71;^ E + 1) = 0(|a|) - 0(ln fc). 

The next equation is 



i4d- 



F P ■ ■ 2d I Afc A2fc A'-i/c, 



which implies 

(60) <,.-J_ft,(?flti).0(|„p 

With this choice of coefficients we get 



1 f h2{a) h^{o:) hr{a) 



2d V k Afc A'-2fc 



= 



fli 



.fc3 fc4 7 ' V ' fc 7 Vfc2 ' fc3 ' --V 2d V Afc A'-2fc 
Considering the coefficients for fc^^ -^^g get 

-a3 + aia2^ 3 — = 03 = 0102 5—. 

a an tt tt^ 
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We get 



2a+l\ 1 , {2a+ 1 



It is clear that in this way we get all Oj, j < r dependent only on {hk}k<j and with this choice 
Xo/k = Tr/d + k~^ai + k~'^a2 + . . . + k~'^ar satisfies 



C. We write h = 1/A. In the general case F''{i;h) = Fq{i) + hFi{i) + ... + h''Fr{i) -. 
Fo{i) + K{i\ h), with Fo(^) as before, Fj polynomial in t of degree (r — j). 

Then F(0; h)/h = K{0; h)/h = i^i(O) + /iF2(0) + . . . + h'-^Fr and we can decompose 

F{hv; h) imh) , , Hjhy) , K{hy;h)-KiO:h) 

—f^=^^+^^y^^'^ + ^^ + h = 

K{Q-h) ^, . H{hy) , . 

where kj is a polynomial of degree < j and fcj(O) = 0. We get 
F (|; 1/a) = ^^(0; 1/A) + -G{y, m) + ^ 

3=2 

Combining this with (|57p where ^j+i is a homogeneous polynomial of degree j + 1, we get 

F (|; 1/a) = F(0; 1/A) + \G[y,t^) + ^ A-^/,(2/), 

where for 2 < j < r, fj{y) = hj{y) + kj^i(y) is a polynomial of degree < j. We have for example 

/2(j/)-Fi'(0)y+iF^'(0)yy. 

We need to solve the equation 



where we write fjia) = fj{{2a + \)/2i) = 0(|ap). 

Then we can get all coefficients in the expansion ^^^ + ^ + ^ + -- - + ^ such that 



\ *=i ^ i=2 



First we get 



1 " 

"1 = "7I;^IZ''^(2a, + l)- 



Fl(0) 



m^"^^ ' 2d 

1=1 
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Let Fj{0) + fj{a) — qj{a) and fix 2 < to < r. Suppose that we have aheady chosen ai, . . . , am-i 
such that 



^ i=l ^ j=2 I 



The left hand side of ([5^ is then equal to 

"2 as , , I 1 92(a) g3(a) 9r(a) 



fc2 ' fc3 k'- ' 2d\ \k \lk A'-ifc 

Denote 

N_a2 , as , , flm-i , 1 { q2{a) , (73(a) ^ ^ g,„_i(Q;) 
Then equation (I62p writes 



7— m 7— m U 



2d ^ Ar'fc V ^™ 



j—7n j—7n 



i^m-i(Ao) T- + O , „ , . as |a|/fc 0. 



Let /(oi, 02, . . . , a„i_i) be defined by 

, . . . , a„t_x) , I 

Then equation (|63p writes 

/(ai,a2,...,ar»-i) g/^ ^^l"^^ ^ 1 gm ^ gj^a) 

fcm ;,m+l y + f.rn + J + 2dXl]'-^k .^^^ \{'^ k' 

We chose such that 

1 / ;\ m — 1 1 / 7\ m—l 



/ + am + ^(^-j (jm(a) = a™ = -/(ai,a2, . . . ,a„_i) - — l^-j gm(a). 

With this choice Aq satisfies equation (|62p with the right hand side O ^ ^^L+i ^ ■ 
We summarize in the following theorem. 

Theorem 5. Lei A^^s be defined in I118\) . For r G Z+ Zei -F^ be defined in |i7|). T/ien /or any 



(fc,a) € N X N" SMc/i that \a\ = ©(Infc) and k large, there exist functions aj = aj{a) — 0{\a\^] 

polynomial in a, j — 0,1, ... ,r, such that if we denote 

(64) 

, / ai 02 ar+l\ TT 1 /r, N -^i(O) 

A. ■.^k[a, + - + - + ... + ^)eA^^s, ao = -, a, = -—^ ^^,i2a. + 1) - . . . 

1=1 

then we have 

2dXr + XrF^ ( ^^ii ; 1/ A,) - 2^^^ = ^ ^ ^^"^ 



2a^ 7 V 

and i/ie coefficients aj, are independent of r for j — 1, . . . ,r + I. 

Note 1. (1) For any r,k CzN large and a — 0{liik), Xq :— Xr{a,k) satisfies 

1 - e-^^'^'-KiXo) = O (^^1^) , n^here Kl{\) e^^^^C^^^). 
(2) Using that operator M is microlocally unitary we know that i^i(O) is real. 
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(3) For large k and a = 0{\nk), we have ReA(a,fc + 1) - ReA(a,fc) = 7r/2 + ©(l/fc). For 
any (3 = Oilnk), j3 ^ a, we have 



1 . ,^f\T?k 



luvXia, k) - IniA(/3, k) = —a ■ {a - [3) + O [ —— , k » I. 

2d \ k J 

In order to have good separation of strings as /c —> oo ( the first term in the above formula 
is dominated over the second error term) we need to impose the Diophantine condition: 

(65) a ^ /3, |/3| < m \fia - > ^(^e"^", D > 0. 
It would imply for any a ^ (3, 

\al m<cl.k^ \,a - > > » O {^^) , cD < 1. 

As the next term is of order k^^ we have separation. 

9. The first local Grushin problem in Wq. 
9.1. Notations. In flint for any G N we can chose r g N such that we have 

Monomials x" are formal eigenfunctions for e^^^^ (-f;!/^). Let Ka denote the corresponding 
eigenvalue: e~'^^^ {i;'^/^)^'^ = Kax". We have 

^ n r 

(66) = exp{-- ^/i,(2a, + 1) - zFi(O) - A'^+i (F,(0) + /,(a))}, 

(67) where fj{a) is polynomial of degree < j. 

We can have Ka{X) — Ka'{X) for a ^ a' . We chose a value Aq G h.A,B and a multi-index ao such 
that 

(68) l-e-^2dAo^^^(^^)^0. 

Let /? G N" be such that KpiXo) ^ Kaoi>^o). 

We suppose that the Diophantine condition ([65]) in Note [T] in the previous section be satisfied 
and write |e-*2''^oii:/3(Ao)| = el''("»-^)l+''. Then if |/3|, |ao| < C In Re A the Diophantine condition 
implies 

(69) |/i(ao -/?)!> c(Re A)"'', r = o ((Re A)"'') , 
with some < (5 < 1. If \^i{ao -(3)\ = ((ReA)-*) then we have 

|1 - e-'2<iAo^^(^^)| ^ ((ReA)-^) , < ^ < 1. 
Let A belong to a small neighborhood SIaq of Aq. 

We chose sufficiently small such that for all (3 G N" with Kp{Xq) ^ Kag{Xo) we have 
A G f^Ao n A^,B, |1 - e-'^^^Kf^iX)] > eo and |l - e-'^^^K^^i^o)] < eo, 

with 

(70) eo = C((RcA)-''), 0<S<1. 
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For A e r^Ao ^ ^A,B we have 

Im A = ^^^{2ao + 1) - ^^^^^ + O {{Re A)-*) , 
where we have IniFi(O) — (see Note[T|). Let 

J = {a e N"; |a| < ClnReAo, Ka,{Xo) = ^^00(^0)} and a = Card J. 
Then we have 

(71) a < 0(l)lnReAo. 
We can chose A^^B(eo) C r^Ao 1^ ^a,b such that 

(72) for A e A^,s(eo), « G J, |l - e-'^''^K^{\)\ < 2tQ. 

We win use the foUowing notations: Ai = ReA, h~ 1/Ai and s = hh\{\/h). 

9.2. Monomials form an almost ort honor mal base in For CC C" we denote 

H^{n) = {u e nol{^) such that / |upe"2*/''L(da;) < 00}, 

Jn 



where $ is the modified weight defined in Section[7]such that $ — |a;p/2 for x € i?(0, e-\/ft,hi(l//i)), 
e < Coo- 
Since the weight function $ is fixed from now on, we shall suppress it from our notations as 
much as possible. Normalized monomials 

(^„(a;) = c^h-'^'^ih-^'^xY, Ca := (^"a!)-^/^ 

form an orthonormal base in iJ$(j(C") for $0 — 1"^- 

Following [5J we show that faix) are orthonormalized in iJ$(il) with an error of order for 
any € N. 

Lemma 8. For all C B{Q, c), where c is some real number, we have 

for iV > 0, which can be taken arbitrary large by choosing coo large enough, and O is uniform in 
Proof: 



{U]^,\x,\'-^) e-'^^^^/'^L{dx) ^U]^, / \x,f-^e-'-^^^'^L{dx,)- 

Jc 

(nj^il^jf "^■) e-l"l'/''L(dx) + / (n^^ilxjf e-'*(")/''L(dx) = 

|x|>e^ft ln(l//i) J{x<£Q, \x\>ey^h\n{l/h)} 

11-12 + h 



and e is such that 0(a;) = \x\'^ /2 for |x| < eyjh\n{\/h). The first term /i is equal to tt"- a\h^°'^'^" 
We need to show that 



7r"a!/il"l+'' 
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ioT h uniformly in a. As $ > c\x\ for some c > 0, it is enough to show the foUowing: 



{x<£C", |2;|>£^/ilii(l/h)} ^ ^ \C/ Jy>ey^cln{l/h) 



<cl 



L{dy) 



where we mean x" = Hj^^x"' and make change of variables y = 
Next we show approximate orthogonality: 

J\x\<e^yh ln(l//i) 
xen, |a;|>e-y/iln(l//i) 

Ii ~ 0, h ^ 0{h^) by previous calculus. | 

9.3. Grushin problem. We denote Vo — tt{Wo), where Wq is the image by kjj of the domain 
defined in ((22l) . 

For e N, we introduce the expansions: 

TMu{x)= 5] (a!ri(9»(0)x"= ^ (a!ri((/ii/29,rz,)(0)(/i-i/2a:r. 

\a\<N \a\<N 

We denote -.^IniTN, := ICerTN, M{N) = dim(F^). 
Suppose that e N satisfy 

(73) max|a| < A^ = 0(l)lnReAo, 



where Aq is fixed as before satisfying ([68]). 
In addition, we need the expansions 

Tju{x) = ^ («!)-! (a»(0)x° = ^(a!)-i((/ii/2a^)«y)(o)(/j-i/22:)«. 

aeJ aeJ 

We denote Fj :— Imrj, Oj := ICerrj. We have a(eo) = dim(F7). 

We have introduced the approximately orthonormal basis of monomials: 

^cix) = Co,h^''^^{h'^^^x)", Ca = (7r"a!)-i/^ 

Then 

ia[)-\{h'^'d,ru){0){h-'/^xr = • inhr/'{a\)-'/^{{h'/^d,ru){0). 
Let i?- : C 1-^ Fj, i?+ : Fj i-> C° be defined by the formulas 

U e Fj, U- = {U-,a}\aeJ £ C", = ^ 

aeJ 

R+u - |(7r/i)"/2(«!)-i/2(/,i/25^)ay(oO ^ 

By direct calculation we get R+u — {u, fa)ae,j{^ + 0[h^)). 
Let 



2 

aeJ 
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(74) 



We write ||.|| for natural operator norm. We have 

R-]-R— = Id^£aj R—R-]- = Tj = Id^pj, T]\fR— = i?— , R^Tpf = 

||i?_|| = l + \\R+\\=l + 0{h^). 

The estimates on the norms of operators R± follow from Lemma [H 
We pose the first Grushin problem for A G Aa.b{^o)'- 

Tn{I - e*2'iAg-«AF-(/;l/A))^^y _^ ^ 

R+u = w+, u,v £ . 

Theorem 6. For a given G F^ x C and all A G AA,B{io), the Grushin problem |y^[ ) 

has unique exact solution (m, u_) G F^ x C, which satisfies 

(75) + \u-\c^ < O ((ReA)^) + k+lc»), < ,5 < 1, 
where 5 is as in Ii69\) . ft2 CC fii CC Vq aKC? |x|c« V^asJ l^aP- 

Proof: We denote h = 1/ReA, A G Aa,b{£o), z = e*^''^. The solution operator for ([74]) is 

where 

a^J,\a\<N 

= {l-zK^)-^(v,ip^)^^{l + 0{h^)), 

a^J,\a\<N 

E"_{v) = ({7Thr/^{al)-'^^h^/^d^rv{0)) = {v, (1 + 0(/i^)) = 

El_iv+) = -{{1 - zK^iX))}\^^jv+, 

where {(1 — zKa(X))}\aeJ is diagonal matrix with entrees (1 — zKa-^), . . . , (1 — zKa^). Some of 
entrees can coincide. Solutions to ([7^ are u — E^{v) + £^*J_(u+), u_ = E^{v) + El_{v+). 
Let CC 5^2 CC f^i CC Vb, s = hh\{\/h). We have the following estimates: 

(76) \\E\v)\\l^<e^^ {^hr{a\)-'\{{h''^d,rvm\^ {l + 0{h''))< 

a^J,\a\<N 

<e^^\\v\\l^{l + 0{h^)), ||i?°(t;+)||^,^ =^K,j2(i^(^(/^iV))^ 

ae J 

as El{v+) = R^v+. 
We get 

(77) ll^^llo, = \\E'{v) + El{v+)\\n, < e^' (l + Oih^)) i\\v\\n, + \v+\c^). 
As E^{v) = R+v, we have already proven that 

mv)\i.<{i+oih''))\\v\\i^. 

By (in]), ^ and jTO]) we have also 

\El_{v+)\c^ <Caeo\v+\c^ < 0(1) InRe A • (Re A)-*|i;+|c<., < S < 1. 
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It implies 

(78) \u^\c^ = W|c= + \El_{v+)\c^ < (1 + + li^+lc"). 

Adding inequalities (TfT]) and ([78]) and using ([70]) we get the estimate. 



As the Grushin problem ([7^ is invertible then it is well posed. Let 
PAr(A) := rAr(/ - e e ^'^ '{X))tn. 

By general consideration it follows that : JCer E^_^ ^ K.er Pm{)^) is a bijection. As Pn{^) is 
a family of Fredholm operators depending holomorphically on A £ ^A.B{^f)) then 

AA,i3(eo) 3 A (Piv(A))"' 

is a meromorphic family of operators. 
We have the formula 

{Pn{\))-' - E^W - EliX) {El+iX))-' E^{X). 
Let 7 be simple loop such that E^^{X) is invertible for A e 7. 

The number of poles n{'y) of Pn{^)^^ inside 7 counted with there multiplicity is given by 

^(7) = ^ iPNwr' = i / (^+ (^-+) ' '^^^ 

From the other side, the number of roots of det + inside 7 is equal to 

-(7) = - / %^^^A ^ ^ ftr\iE\y{E\y^}dX. 
As in [8] we have then 71(7) = 771(7). 

10. Some useful estimates. 

10.1. Estimate on un G ICerriy. Let s = h\n{l/h), h = 1/Ai, Ai = Re A, A G Aa.b, B,. = 
B(0,r) = {x, \x\ < r}. 

Lemma 9. Let N eN, N = ©(In Re A), satisfying (W. Let un e H^{B^r^) n Kerr^. Put 



zvsAi 

Let e > rm be such that <f> = |a;p/2 for x G B^^. Then for any tq < rm, 

Proof: We need to prove 

e-l"l'/''|uN(x)pL(da;) < / e'^^^' /''luNix)^ L{dx). 



We change the variables x = i^/i in order to have the domain of integration independent of h : 

X e ^ X e B^. 

Denote $s = ^^{y/sx) and 

H^^{B^) ^nol{B,)nL^^{B,) ^ {v eHol{B,), [ \v\^e-''^'^^/^L{dx) < 00}. 

JB, 
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a! 

\a\<N 

Suppose that e > is small enough such that <& = |a;p/2 on the domain of integration. Thus 
= |a;|V2- We must prove that for {s\i)~^^^ N + 1/2 < e and = l/2(sAi)"^/ + 1/2 
we have 

The monomials constitute an orthogonal base in _ff$ (Be) and it is enough to show the inequality 
for x", \a\ > N. Let aj e N be such that aj > N. For un{x) = x"^ we have 

'|x|<e J r<€\/sXi 

As 

we get that the function e^'' r^"^^ is increasing for < r < r™. Then 

For all a such that jal > N the maximum r,„(a) > rm(a7v) with jaArj = A^. Thus, if the above 
estimate is valid for |aAr| = N then it implies that it is valid for any |a| > iV. ■ 

10.2. Model estimate. 

Lemma 10. As before we write s = hhi{l/h), h — l/Ai, Ai = Re A. Let be any X-independent 
neighborhood of and rg > such that B^^^ CC il. 

Let A = diag(e''i, . . . ,e^") be as in Theorem\^ Denote {A^^x, x £ 9.}. 

For any cq > and \z\ < e"^"'"-^! and ro large enough there exists m < 1/2 such that 

g-2*(x)//.|^g-E?.iw(-.a.+i/2)y|2^(^^) < / e~^^^^'>^''\u\^Lidx). 
Proof: We have 

e-2:"=iP,(^.9.+i/2)y ^ e-^l=it'o/^u{A-^x), A = diag(e^i . . .e'""). 

Let first $ = \x\'^/2. We put for simplicity n = 1, /ij = Then e^^'^^+^/^^u = e^^/2u(e^^a;). 
The general case is straightforward. We change the variables x — e~^x, L{dx) = e^^L{dx), 
and denote the new variable x again by x. Let A~^{n) := {e~^x, x G fl}. Then we have, with 
s = /ilnAi, 

e-f"- f e-l"l'/''|zM(e-^x)|2L(cix) < e2=°'"^i+^ / e-^'"\''\'/''\u{x)\^L{dx) = 

(79) 

<g2colnAi+Mg-(e^''-l)e---r^//. ^ e"l"l'/'Xa;)pL(dx) = 

= m^ [ e-l=^l'/''|w(x)|2L(dx) 
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with 



^2 _ g2coln Ai+Ai-lnAi(l-e-2*')r^ 



We have m < 1/2 if 

g2coinAi+M-inAi(i-e-='')r= < ^ 2co hi Ai + /i - In Ai (1 - e-2^)r2 < -2 111 2 4:^ 



/2colnAi +/i + 21n2 



InAi • (1 - e-2A') ■ 

For general weights $ we use the strict convexity and get (after change of variables x = e~^x) 

$(e^x) - $(5) > c|ip, c > 0, xfzVi, 
and in the above formula (|79|) exchange e^^ — 1 by c. | 

10.3. General estimate. In this section we will write h = l/X and suppose that A is real. 
Generalization to A G Aa.b is straightforward. 
From ([53 it follows that 



where the polynomial Pj{y) = X]|a|e[j+i 2j] lajV'^ is a linear combination of monomials of degree 
in [j + 1, 2j]. Moreover, pi{y) = -^H"{0)y ■ y. Then 

e-'^^'^y^/'^ = e-'^^y (1 + hp^{y) + h^p^iy) + ...) = (l + hpi{td^) + h^p^iid^) + . . .) e'^^^y . 

Taking y — • I — {xdx + 1/2) /z, we have the following representation: 

71 

i=\ 

polynomial of degree r independent of h. 
Then, 

e-^^''^^'^>^)l^e-'>:^'^.=ii'^'^''^^^'-^l'^)u ^ e-*^(*^^)/''e-^^=i^^/2M(A-ia;), A = diag(e^i . . .e^"). 
For simplicity we put n — 1. We have 

r-1 

^e-''''^u{e~''x) + h^h^-^ qi^jdl{e-''/\{e-''x)) = 

i=i ie[i+i,2j] 

= e-^^^uie-^x) + e'f^^T^ie-^x), 

r~l I 

where rj'u(a;) /i^/i^^i ^ ^ Cfc,,M(''')(a;) • a;'= = ©(/i). 

In general case 

^-^F(hy■h)/h ^ ^^^K{0■,h)/h^-^Giy.^.)^^ ^ ^ ^^2^^^^) ^ ^ ^ )^ 

where (/^(y) is a polynomial of degree at most 2j, and as before, this can be written 
We will write it in the form 

F{hy-h) F{p^ , , J{hy;h) 
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and then (see also [13], page 347) 

^~iF(I(h)-h)/h ^ g-j_F(0:/i)//ig-a(j/i9j,:/i)//ig-2:r=iA'i(2^.ft + l/2) 

and the difference from the case F ~ G{I.ji) + H{I) is only the factor g"*^'"-'')/'', which is 
bounded as /i ^ 0, and that qj{y) is a polynomial of degree at most 2j and not just in [j + 1, 
We have 

^~iF(I(h)-h)/h^ ^ ^-^F(0■,h)/h^-lJ(^h^^■,h)/h^-Y.7=l^^ii^i9^ + '^/^2■)^ = 



-iF{0;h)/h \ 



r-1 2] 
]=1 1=0 



(80) = e~'^io-M)/h (^e-^/2u(e-f x) + ^-''/^^''(e-^x)) , 

r-l 2j i 

where T^u{x) :— 



We put z = e~'^^'^^. Then for A G Aa.b we have jz] < e^''^^^^^^, with A as in the definition 
of Aa,b- 

We can apply Lemma [TOl with cq = 2dA. Then, for tq large enough we have 



g-2*(x)//i|^|2 ^_ r''u(e'''x)| i(dx)) < 



(81) / e-2*(-)/''|u(a;)+r'-u(x)|'i(dx). 

JA-i(n)\A-i(S,„yi) 

For any CC ili CC fi, we have for any u £ i/$(ri), 

see Corollary 4.2 in [8]. 

If X € il C Vb then |x| < 0{l)h\n{l/h). Taking u e /Cer(rjv) and using Lemma[9]we get 



\\{h-^/^xnh'/^D,'fu\\n, < O ((ln(l/;i))l"l+l'3|) \\u\\n^^B. 



1-0 v^' 



Thus the partial differential operator is bounded of norm 0(h) : iJ$(ili) i-^ _ff$(f72) and as 
u G ICerirpf) we have, 

(82) llT'-^lk^, < Oih)\\u\\n, < Oih)\\u\\n,\B^^^- 

Then, using ((821) with Qi = Q and = ^"^(Jl), we get 

f e~^^^^'^/''\T'-u{xf L{dx)<0{h') ( e-2*(")/''|u(x)|'L(dx). 

We can choose tq large enough such that 

the right hand side of (EID < i / e-^'^''''^/''\u{x)\^L{dx). 

4 Jn\B _„ r- 

This implies the following lemma: 
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Lemma 11. Suppose A G h.A,B o,nd |A| > C large enough. Let h = 1/ReA, tq, e be as in Lemma 
\^ Let r2 be any X-independent neighborhood of and 

B,^^ CC B^^ CC n. 

Let ^ = diag(e'"i,...,e'"") be as in Theorem^ Let A-^{Q) := {A-^x, xen}. 
If rQ and e are sufficiently large then for some m < 1/2 

(83) / e-2*/''|e-*2''^e-^^^''(^'i/^)u|2L(da;) < to2 /■ e-^''/''\u\^Lidx) 
for any u G H(i,{VL) n /Cer(Tjv). 

10.4. Bounds useful for the second local Grushin problem, Theorem [7l We suppose 
A e Ka,b- Vi fix g N and let un S i?$(Vb) n K,er{TN)- Let ro > be large enough. In all 
lemmas below we use the following convention for the domains: 

(84) B,^^ CC CC ^2 CC r!i CC Fo. 

Lemma |6] implies that the norm of Rr — Mq — e~^^^ is small in H^{Vo). We have 

Lemma 12. For any N £ N there is r — r(N, N, tq) G N sufficiently large such that 



< 



< 



with domains satisfying ^8^. 
Proof: Let z = e^'^rfA^ -y^g 

use Lemma |6] (with N instead of N) and Lemma [9j 



Combining Lemma [TT] and Lemma [T^] we get 

Lemma 13. Let Mq = e-'^-^'t-^-i/^) + i?''. Suppose A G Aa,b, a«rf ^-^(r^) = {A-^x, x G 
fi}, A = diag(e''i, . . . ,e^"). 

for any N Cz N there is r G N swc/i that if ro is sufficiently large then for some ni < 1/2 

(85) l|e-*''^'Mo^A'll^,\s,^^ < m2||u^^||2^^^_,(5__^^) + Or{h^)\\uM\l,\B^^^ 
and 

(86) ll"^llo.\i3,.„^ < Oiimi - e-^^''Hlo)uN\\l,\B^^^ + Or{h^)\\uM\\l,\B^^^. 



Proof: We show (|86l) . Choose e such that <i> = |2;p/2 for a; G B^^ and tq < 
before. Then, with z ~ g-«2dA^ ^^^^ some constant to < 1/2, bound [85l implies 



2VsAr 



as 



> 



> \WN\\n2\B,^^ ~ "1 \\uN\\n^\A-HB,^^) - ^r{h )\\uN\\n^\B^ 



(87) = (1 - m'^)\\uN\\l^\B^^^ - m^\\uN\\l^^^\A-^B,^^) ~ Or{h^)\\uN\\l^\B^^^- 
Lemma ^ implies 
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Then we get that 

the right hand side of (Ell) > (1 - 2m^)\\uN\\l,\B ^ ' Orih^)\\uN\\L\B ^- 
Then we have 

,2ml,_j|2 ^wij „T,t\„._A\2 I /f.JVMi j|2 



Lemma |9] imphes that |luAr||^ < HujvUo , „ which together with estimate ((86)) imphes 

Lemma 14. Let A G A^.s- For any G N i/iere is r e N such that if ro is sufficiently large 
then 



(88) h^ll^, <0(l)||(/-e-*^'''^MoW||^,+a(/i^^)hA'||^,\s,^ 

wif/i domains satisfying [8^. 



(89) 



11. The second local Grushin problem. 

With the same notations as in Mq = e-»^-F" (J:i/A) ^jir^^^ Aa,b{£o) verifying ([72]), for 
any Q2 CC CC Vq, we pose the second Grushin problem: 

(/ - e~''^'^^Mo)u + R^u^ ^ V, 

R+u w+, u e H^{VLi), V e H^{^2), u^,v+ e C°. 

Theorem 7. Let A G A^.B(eo). We consider the Grushin problem /189\). 

Then for any N Cz N there zs r G N large enough such that we have 
(90) 

WuWn^in,) + |«-|c» < O ((ReA)^) (Mn^^n,) + \v+\c^) + Or (|RcAr^) \\u\\H,(n,^, 0<S<1. 
Here S is as in \69fl . 

Proof: Let z = e"*^'^'*', h = 1/ReA, and A^ sufficiently large. Applying (1 — tn) to the first 
equation in (|89l) and using (1 — tn)R- = 0, we get (1 — tn){I — zMo)u = (1 — tn)v which implies 

(91) (/ - zA/o)(l - tn)u = (1 - tn)v + tn{I ~ zMo){l - tn)u - (1 - rjv)(/ - zMo)tnu. 
Then we apply estimate (|HH1) (LcmmalH]). Let B^^^ CC VI2 CC Oi CC Vq. Then 

11(1 - ta,)u||o2 < c||(l - tn)v\\q,2 + 

(92) + 0{1)\\tn{L - zAfo)(f - tn)u - (1 - rA,)(/ - zMQ)TNu\\n,+ 

+ a(/i'^'')||(f-rAr)«||oAB,„^. 

We need to bound ([92| . 
We have 

||rArze-^^^''(^^i/^'(l - ^Jv)"llo. < 0(/i^)hlb,, 

where CC CC 51i CC Vq, and similar, interchanging t^t and 1 — t^t. It follows from the 
following facts: 

if u G Oiv = /CerTA,, then m^'^) (e-''a;)(e"^a:)^ G Oat V/3, 
and formula ([SO)) , which implies that Tjvze^*'^^'^(^'^/'^)(1 — rjv) = for any fixed order r. 

Then it is enough to estimate \\t^zR{1 — tjv)||o2 ^-nd II (1 ^ '''n)zRtm\\q,2 which can be done 
directly, but from more general result, Lemma [6] on Vq, we know that for any N we can find r 
such that (HD) < Or{h^)\\u\\^ . 



SCATTERING POLES NEAR THE REAL AXIS FOR TWO STRICTLY CONVEX OBSTACLES. 



39 



We get the estimate 
3) 

Apply TN to ([89]). Then 



(93) \\{l-TN)u\\n,<0{l)\\{l-TN)v\\n,+Or{h'')\\u\ 



tn{I — zMo)u + R-u^ = TpfV 
R+TiqU = v+ 



since tjsiR- — R- and i?+ — R+t^. Then 

R+t^u = v+. 

The first term in —tn{I — zMo){l — tn)u — tnzR^tnu which is equal to rjvzMo(l — tisi)u = 
TNze~^^^'^^^'^^^\l — tn)u + tnzR^{1 — rjv)w is already estimated. The second term can be 
bounded as before using Lemma [6] in Section [6l 
Then estimate (|75|) implies 



(94) \\rNu\\n, + |u-|c« < O {h-') i\\rNv\\n, + \v+\c^) + Orih^)\\u\\^^. 
Adding ^ and ^ we get 

(95) \\u\\n, + |w-|c" < 0(1)I1(1 - mMn., + O {h-') i\\mv\\n., + k+|c») + Or{h^)\\u\\^^. 
Then using Proposition 4.3 in :8], we get the desired estimate (|M| . | 

12. The global Grushin problem 

We denote fiint G T*dfli the A-independent neighborhood of ai such that Wq CC ilint CC Wi, 
where Wq, Wi are defined in ([22]) . ((26)) . Let ficxt be such that Cficxt CC fiint- 
We consider the original operator / — M(A) in Ti. equipped with the norm 

\u\h ■■= !l(l - X2)TA,Gw||L2^(T-(af2i)) + llxit^u|lL2(e-2*/'>L(da;)) =: l"|n„,t + |w|n.„t> uen, 

where xi G C'o" (*C^") is equal to 1 in a neighborhood of ai = 0, with the support independent of A 
and Wo CC suppxi CC riint CC Wi, and X2 G Co{T*X) equal to 1 near (ai,0) and essentially 
the same function as xi after suitable identification of domains. Here L^q{T* (dili)) stands for 
L2(A;e-2^/'^|(a4)|2™da). Let 

U: ff*(l^i„t), U = S.FTBargman, Vi D neigh(ai) ^ neigh(O) C Vi, 

7r(Wi) = Vi, be as before (we omit') such that 

UM{X)u = e-^"^^MoUu + 0{h'^) \\u\\n, in -ff<i,(f^int), Mq - e"'^-^" = Oip'^''+\ h'^''+^), 

where h = 1/ReA, A e A^^^. 

Let 1^ : i7$(ili„t) i-^ -ff (Ate) be the approximate microlocal inverse of U such that if x(a;, hD^) 
is pseudodifferential operator adapted to H{AtG) with compact symbol and with suppx CC r^inti 
then xiVU - /) and [uV - I)x are neglectible, VUu -u = 0{h^)\\u\\-H. 



Theorem 8. For a given (u, w_|_) G Ti. x C" consider the Grushin problem 
(96) 



(/-A/)M + Fi?_U_ =w, 

R^Uu = 



40 



ALEXEI lANTCHENKO 



For X €z A^^B(eo), verifying |7ii[ j, and Re A sufficiently large the Grushin problem 1196]) has a 
unique solution (u, w_) G 7i x C°. Moreover, we have the a priori estimate 

(97) \u\n + k-lc" < OiiReX)') {\v\n + \v+\c') 
with some S, < S < I. 

Proof: The existence of a solution follows as in jOj . We need to show estimate (|97p . We apply 
U from the left in the first equation in and denote Uu — u. Then 

r (/ - e-^"^^Mo)il + R^u- =Uv + w, w -.^ (UM - e-^^'^^MoU)u + {I - UV)R-u_ , 
[ i?+{t = V+. 

Denote v = Uv + w. Introduce as before fJoxt, f^int, Wo, Vb — 7r(Wo). Let 

CC r^i CC Vb CC CC CC VL\ CC f^int- 

Estimate on Vq : 
By dMl) we have: 

(98) \\u\\H,(n,) + \u-\c^ < O {h-') iMH.in,) + k+|c») + a {h'') Mn^in,), < S < 1. 
We have 

(99) MH,(n,_) < Or{h^) {\\u\\H,in,) + |^^-|c») ■ 
This gives 

(100) Mh^hw + \u-\c^ < O {h-') {\\Uv\\ 
Estimate on flint \ Vo ■ 

In (/ — zMo)u = —R^u^ +v we use (a variant of) Theorem[3l estimate (|50p . which implies that 
there is A^i > such that 

lkx^^ou||L|(nj,) < 0{h^^)\\u\\H^(n{), suppx C r^mt \ Vq. 

We have then 

||x(l - zMo)m||l2 > ||xw||L|(oi) - 0{h^^)\\u\\H^(n{) 

and 

llxwIliKni) < \\xR-U-\\Ll{n'^) + +C'(/i^')ll"llff^.(n'i)- 
This implies for any g N 

(101) IIx"IIl|(o;) < 0{h'')\u4c^ + llx^llLl(ni) + 0(/i^^)||u|k,(ni)- 
Similar, applying (jSOp to each term in w separately, we get 

(102) WxwWLliii',) < 0{h''^) (lIx^Llaii) + \u-\c^) 
and we have 

(103) llx^i||L|(o;) < 0{h''^)\u.\c^ + WxUvhiin',) + 0{h^^)\\u\\H,in[). 
Applying (I1GG|) and ()103|) with u replaced by Uu we get 

(104) \u\n'^ + |u.|co < 0{h-') {\v\n'^ + \v+\c^) + O (/i^^) {\u\n'^ + |u.|cO , 

where fij CC ft'i CC i^int- 
Estimate on f^cxt : 

In order to estimate u — Mu — VR-u^ + u we use bound (I42p with some > : 
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and get for any N E N 
Summing up we get 

|u|f2„t + |w|f2i„t + \U-\C- < 

with A^o ~ min(iV, A'^i, A'2). The last term can be absorbed in the left hand side and we get 
estimate (|97l) . 

Then it is clear that 

): TixC^^HxC" 

is Fredholm of index 0. Estimate implies that V is injective and thus bijective for A € 
AA,s(eo)- 

I 

Denote 

E E+ 



^ ^ ' E^ 

the inverse of V. Let h — 1/Rc\, X E A.a,b- Then it is known 

e cr(/- Af(A)) e cr(S+_(A)). 

Let v+ e C°. If iu,u-) = £(0,u+) then u = E+{v+), m_ = E+_{v+). 
Let 

r {I - M)u + VR-u^ =0, 

Then we have 

(/ - e-2*''^Mo)u + w = ([/Af - e-^''^^MoU)u + (/ - 

and w satisfies 

< ai/i"^) + |u_ic") = a(/i'^)k+ic«. 

We have then for r large enough 

Tj {I ~ e-2''iAe-*^^''(^(i/^)'i/^)) Tju + R^u^ = Or[h'')\v+\c^, 

TjU — 

and 

(/ - e-2'<iAg-^AF'■(/(l/A);lA)^ ^ ^ Orih^)\v+\c^ ^ 

= -R+ (/ - e-2»'^^e-*^^"(^(i/^)^i/^)) R-V+ + Orih^)\v+\c^. 

We have then for A G Ayii.s(eo) : 

i?+_(A) = + a ((ReA^^,)-^) , where EI_ = - {{l - e^''''K^{\))} ^^^^ . 

Let 7 = 7(AJ^ be simple loop such that i?1_|_(A) is invertible for A G 7. Suppose 

Ae7 A = A^+a((ReA^^,)-^), 
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where f. is solution of ([55]) . Then the number of poles of (/ — M{X)) ^ inside 7 counted with 
multiplicity is equal to the number of roots 77^(7) = '7i(A^ j.) inside 7 of det E^j^(X) — 0, 

(det£;0+(A))' 



Ml) = 77- [ 



det 



-dX. 



We have G - M(A)) if A = A^ + a ((ReA^ fc)"^ 

This accomplishes the proof of Theorem [21 | 

Appendix A. Some facts about the quantum billiard operator M. 

A.l. Boundness. Let fti be an obstacle with analytic boundary, non-trapping. Let J7i C 
B{0,R). Denote 

(105) Ua,b ■■= {A; ImA < A\X\^/^ - B}. 

Let R{X) be the outgoing Dirichlet resolvent in Cili. 

Theorem 9. Let x G C^f (M") &e such that x = 1 m B(0,i?). T/iere esisf A, S > smc/i f/iai 
Xi?x('^): defined for Im X < /las an analytic extension to the domain U05]) as a bounded operator 
L^(Crii) i-^- -ffQ(Cr2i), satisfying the estimate 

3C,D> 0; ||xi?xll£(L=(Cn,),i/,5(to,)) < 
where Im A"*" = max(Im A, 0). 
For i = 1,2 let 



be the outgoing resolvent of the problem 

(A + X'^)Hi.+ {X)u = in M"+i \ Qi 
Hi^+{X)u\Qi2^ = u 

extended as an operator H^^^{dn,) ^^ Hl^^iCQ,). In T it is proven that this resolvent, analytical 
for Im A < has an analytical extension to the domain of the form (jlOSp as a bounded operator 
H^^^(dVli) ^ _ffQ [Qj,(Cf2i) and satisfy the following estimate: 

Vi? > 0, 3C> 0, 3D > 0, VA e Ua,b, mAc(mnion,),miZa,nBio,R))) < Ce'''^^\ 

We denote H2i{X) — i?2,+ o 72 ° -^1,+, where 7^ is the operator of restriction to dfti, the 
outgoing resolvent of the problem 

(A + X^)H2iu = 0, in CrJs 

We define Hi{X)u = Hi_+{X)u\QQ_.^^, where d^ly, — d^li and 

Af(A) = i/2(A)ffi(A) = 7ii72,+72ifi,+ . 

15 (Burq). Operator M{X) defined on H'^{dQ,i) ^ H'^{dni) for ImA < has an ana- 
lytic extension in the domain of the form I1105\) and there satisfies the following estimate 

3D >0,3C> 0, VA e Ua^b, ||M(A)||£(Hi/=(ani)) < C\X\'e^'^^'' 

Moreover, as in a neighborhood of fli, Mu satisfies (A + }?)Mu — 0, we get 

3D > 0, Vs e K, 3a > 0, VA e Ua,b. \\M{X)\\c(H^(Qn,)) < Cse""'^'^^ \X['+\ 
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A. 2. Microlocal unitarity of M with respect to the flux norm. We suppose that A is real. 
Denote h = 1/A. Let v satisfy 

, . r Pu (-/i^A - l)u = 0, u outgoing 

[ 7iii — V. 

We define : i— > I?'(Cf2i) the solution operator. In the similar way we define 

H2^+ : V'{dD.2) i-> X''(CrJ2), which satisfies PH2,+v = and J2H2. — V. Let Hi = 72^?!,+ and 
H2=-fiH2,+ .Then M ^ H2H1 : H^{dni)^ H^{dni). 

The billiard operator M can be identified with the monodromy operator M as in [24] in the 
form presented in ([H], p. 360). Let keTp^{P) be local kernel P near pi = (ai,0). The solution 
operator K in ([II]) is the operator 

e kerpj(P), i/i^+uiaf^^ = u. 

We have 

M = 7Ji,+7iiJ2,+72 : kerpi(P) kerp,(P). 
As in ([S]) we have identification kerpj(P) ~ 'D'{dili) via iiT = . Then the monodromy 
operator M on 2?'(9r2i) satisfies 

ifMw ^ MKv, V e x>'(arji). 

Let X G C°°(Cf2i) be a microlocal cut-off function supported in a neighborhood W2 of fii such 
that X = 1 in Cf)i n W^i, where Wi CC VF2. 

We define the quantum flux norm (see [14], p. 360) on the outgoing solutions u of (|106p as 
follows 

IhllQF (;^[P,1-XM«). 

It is easy to see that |1u||qf is independent of x which implies that AI is microlocally unitary 
with respect to ||.||qf for real h : 

Lemma 16. The billiard operator M is microlocally unitary for all real A with respect to ||.|1qf : 
IjA^i^wHqF — ||i?i,+Aft;||QF = \\Hi,+v\\qf + 0{h°°), where u — is the outgoing solution of 

ilOb]) and WFh{v) C neigh(/3i). Here the wave front set WFh{v) is defined as in [9\. 

Using the Green's formula in a neighborhood of ili bounded by dili on one side, we get 

II«IIqf = ( 1 - X]u\u) = ^{{h^A + l)xu - x{h^^ + l)u\u) = 
h f 

= - / {{d,Hi,+v)W^ - Hi,+v{d,W^)}Sidx) = 
* Jani 

= 2Re {^d^Hi^+v\v) L2(^g^^^ = 2Re {Av\v) l2 (^qq^^ , u^Hi,+v, 

where d^, is the normal derivative, Av :— jijd^Hi^^v. As A is pseudodifferential operator of 
order 0, elliptic near (ai,0) (we use the explicit WKB construction of in the hyperbolic 

zone) then we have for v with WFh{v) C neigh(pi) 

\\Hi^+v\\qf ^ ||w||L2(anjnnoigh(pi))- 

Lemma 17. There exists pseudodifferential operator B of order elliptic near ai such that 
BMB^^ is microlocally unitary for all real A with respect to \\-\\L^{dQ,i) ■ 

\\BMB~\\\L2ion,) = WF{v) C neigh(ai,0). 
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